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ABSTRACT, 


A new general method of calculating the high resolution 
NMR spectra of all molecules with twofold frame symmetry 
has been developed. Full use is made of all of the 
symmetry factoring of the Hamiltonian matrix resulting 
from the twofold frame symmetry, and from any magnetic 


equivalence. 


Explicit formulas have been derived to allow the 
direct calculation, from spin quantum numbers, of all 
of the Hamiltonian matrix elements. The formulas are 
valid for molecules with or without magnetically equivalent 
groups of nuclei, including the important class of molecules 
containing symmetrically equivalent pairs of magnetically 
equivalent groups. For molecules without frame symmetry, 
the formulas reduce to those of the Composite Particle 
method. They are suitable both for hand calculations 
on small molecules and for computer calculations on 


large molecules. 


A systematic three stage approach to the factoring 
of the Hamiltonian matrix has been described, in which 
this symmetry factoring, molecular total spin factoring, 
and species total spin factoring are applied in succession. 
It has been shown that this approach to factoring simplifies 
the derivation of transition selection rules, and leads 


directly to a description of the structure of an NMR 
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spectrum in terms of subspectra. 

A computer program, NUMAR, has been written to 
allow the application of the full three stage factoring 
scheme to the spectra of molecules too large for hand 
calculations. In addition to the usual energy level 
and transition lists, the program provides all of the 
quantum numbers necessary to define the structure of 
the calculated spectrum in terms of subspectra. The 
efficient design of the program permits calculations 
on many large symmetric molecules whose spectra previously 
could be calculated only by approximate methods, if at 


all. 


The new factoring procedure has been applied in 
the analysis of the NMR spectra of several molecules 
with twofold frame symmetry, including the perfluoro- 
2-butenes, the 2,5-dimethyl-3-hexenes, and [(F,C),P],NH. 
In each case, initial NMR parameters obtained from a 
partial or approximate analysis were refined with the 
aid of exact spectra calculated with the NUMAR computer 
program. The use of all of the factoring appropriate 
to each molecule simplified all the analyses; in the 
case of the 2,5-dimethyl-3-hexenes, it permitted accurate 
calculations on spin systems too large for previous 


computer programs. 
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CHAPTER i] INTRODUCTION 
ES? High Resolution NMR Spectra 


In a Nuclear Magnetic Resonance experiment, a 
sample containing magnetic nuclei is placed in a strong, 
uniform polarizing magnetic field and a weak radiofrequency 
magnetic field is applied. For particular frequencies 
of the radiofrequency field, resonant absorption of energy 
occurs; the intensity of absorption as a function of 
frequency is «the "NMR spectrum ‘of the sample, ‘This ‘Thesis 
is ‘restricted to High Resolution NMR, where the sample 
tsea mobile, isotropic fluid, and there are no strong 
couplings to paramagnetic or quadrupolar species. The 
Soeourum then normally consists of sharp,ylines with a 


Maio ei anewioth ot Lt -iHz cor eae. 


Most High Resolution NMR spectra are obtained 
PROMEDUGLeT OL spin. = 4. .- Because nuclei. of .spimn | =.0 


(such,.as 47C and 1°0). have no magnetic moments; they 


t For general information on the nature of Nuclear Magnetic 


Resonance, see Pople, Schneider, and Bernstein (pop. 

For a more mathematical treatment of the theoretical 
aspects of NMR, with proofs of many useful theorems, 

see Corio” especially Chapters 5 and 2. The notation 

Poe imitmein ioe esis wecteextension of Gorio's, is. described 


in Appendix A. 
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do not interact with magnetic fields, and have no NMR 
SUCH owe euucl lei OLospint Los 2 -are normally coupled by 
their electric quadrupole moments to fluctuations in the 
Pecal clectiiicericld cradients, resulting from molecular 
motions. As a result, the lifetimes of their spin states 
are short, and their NMR spectra have linewidths orders 
of magnitude too great to be considered high resolution. 
When a nucleus of spin % is strongly coupled to a quadrupolar 
nucleus such as ‘*N, its lines may also be broadened 

(see PSP. More often, nuclei of spin greater, than % 
may be treated as non-magnetic, because they have no 


Observaple etfect on the spectra of the spin % nuclei. 


Nuclei with electric quadrupole moments are not 
excluded completely. from this Thesis, because they may 
yield High Resolution NMR spectra when they are in 
environments of high symmetry. However, the main emphasis 


Peeounucte, ofsspin %. 


The parameters of chemical interest which determine 
tiestuseh kResolutzon NMR-spectrum of a molecule are one 
chemical shift V5 for each magnetic nucleus’ and one 
coupling constant wee Connecting eacin pair Ol NUCret. 
Porenenicicwuthere are thus a total of n(nt1)/2 parameters. 
The chemical shifts vary linearly with the strength of 
Pieeo st erie ine Magnetic t16eld but the coupling constants 


are independent of field strength. 


There is no practical general method of calculating 
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the parameters directly from the observed spectrum. Instead, 
one must calculate the spectrum that corresponds to a set of 
trial parameters, and compare it with the observed one. 
If the observed and calculated spectra agree, the trial 


parameters are probably correct.! 


in principle, the calculation of a High Resolution 
NMR spectrum 1s quite easy. The Hamiltonian’ has the 


Simple form 
Hise}y, 0) (2) ern] es, Sh) -iciy. (1:1) 
1 LS) 


A suitable basis set consists of products of spin functions 


sorrthesindividual nucleiynoftthe form 
M |I,,m,>. . (1. 2a) 


t ties alWays pOssible=to invert. the signs of all of 


the coupling constants simultaneously without in any way 
affecting the appearance of the calculated spectrum (see 
Corio’, Delos) Ortenetie, spectrum 1S ealso insensitive 
to the relative signs of many of the coupling constants. 
One method of determining the relative signs of coupling 
Sonstantssis described in §3.3.2, The absolute signs 

of some coupling constants have been determined from 
spectra of partially oriented molecules in anisotropic 


media.> 


Throughout this Thesis only nuclear spin Hamiltonians 


are used. The nuclear spin operators such as ill and i are 
dimensionless. See Appendix A for definitions of all 


operators, eigenvalues, and basis kets. 
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to -I, there are 
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basis erunctaons in the set. 


A matrix representation of the Hamiltonian (1.1) 
Mmuewioebasts (liza) is: constructed: ats: dimension is 
given by (1.2b). When the Hamiltonian matrix has been 
diagonalized, differences between its eigenvalues give 
the transition frequencies. The corresponding relative 
intensities are given by the absolute squares of the 


matrix elements of the transition operator’ 
2 Staph oe Bly fr (de 3} 
1 


Pieemecetx ehementsyor X* are most?easily*calculated?in 
Evemeinpre=proauct= basis’ (1.2).>° The matrix elements 

of X in the basis in which H is diagonal can be obtained 
from them using the unitary transformation that relates 


Preevcscs.  “The“total~spectral “intensity that is 


t This equation is an approximation which introduces 
an error far smaller than the usual uncertainty in 


expemimentabeintensities. See Corio” Di oa 
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obtained using (1.3) is? 


wd 


{ ae aitche Papert eae (1. 4a) 
ah 


which for n nuclei of spin % reduces to the simple form 


awa (1. 4b) 


The method of calculation described above is valid 
Potpedds Spantsysitems,, bute practical only for small-ones. 
Porecxomple. seach of the-perfluoro-Z-butenes (55.1) has 
erent mucle1 of spin. %.. The 256x256 Hamiltonian matrix 
has. 32896 distinct matrix. elements. (many. of.which are, 
Diet, 2670). Ihere could be as many as 11440 transitions. 
Pmt Nousanase Ob distinct, transitions;were visible an 
the experimental-spectrum, there would be little hope 
Brerecoonizineg patterns helpful: in estimating the parameters. 
Pr aeeond error, calculations, withithirty-six independent 
parameters, would also be unlikely to yield a correct 
Solution, even if reasonable values for some of the 
parameters could be estimated from previous work with 


related compounds. 


The observed spectrum of each of the perfluoro-2- 
butenes contains only about 100 distinct lines of Silent tcant 
intensity, largely as a result of the high symmetry of 


the molecules. 


For each isomer, the six -CF, fluorines have one 
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chemical shift, and the two =C-F fluorines have a second 
chemical shift. Each spectrum consists of two multiplets, 

one centred on each chemical shift. The relative chemical 
shifts can be varied by using a variety of polarizing 

mierda strengths.**Forcsypical field strengths! {overs 10 kgauss), 
the two multiplets are far enough apart that the appearance 

of each multiplet is independent of the relative chemical 
shifts. Thus none of the eight chemical shift parameters 

need be determined in order to calculate the pattern 


of linesyan,.the multiplets. 


The symmetry of the molecule also reduces the 
number of independent coupling constant parameters from 
Zor toe Pour (F,C-C=C-CF 


BRC=CUR}> FC (2) CRA PPHC=C> OFL)e 


3? 3°? 
Thus only four parameters, rather than 36, need be found 
in order to completely analyze the spectrum. Quite generally, 
the complexity of an NMR spectrum is more closely related 


to the number of independent parameters than to the 


number of nuclei. 


The major goal of this@thesis 1s to outline a 
method of calculating NMR spectra, such that the same 
factors that simplify the observed spectrum also simplify 


thneecarculation. 


t As each multiplet is symmetrical about its chemical 


shift, the two chemical shifts can in fact be obtained 


by inspection. 
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A ewe Factoring the Hamiltonian Matrix 


The NMR spectra of the perfluoro-2-butenes, and 
other molecules, are simplified by two molecular properties: 
symmetry, and large chemical shift separations. As a 
result of these molecular properties, certain symmetry 
operators and spin operators commute with the Hamiltonian. 
This in turn leads to the simplified method of calculation 


thateis the goal of this Thesis: 


If two operators commute, there exists a set of 
vectors which are simultaneous eigenvectors of both operators. 
Piteempiles sthat 1f°a basis) tor<the Hamiltonian has been 
Colmetructed Out of cigenfunctions of any operator O which 
commutes with the Hamiltonian, there can be no off-diagonal 
matrix elements of the Hamiltonian between basis functions 
Delonging to different eigenvalues of 0, ‘The complete 
Hama.lionilan matrix, before diagonalization:, can be reduced 
to block diagonal form, with each block corresponding 
omamcis ferent eigenvalue of;0. If there is an additional 
operator that commutes with both H and O, these blocks 
(submatrices) can themselves be factored into blocks 
eionemthe main diagonal, and so on. Each block is labelled 
by a unique set of eigenvalues of the operators that 


commute with Ht. 


Tian NNinecalculation,. this factoring process ‘can 
be viewed as taking place in three stages, each stage 


being associated with a different kind of operator that 
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commutes with H and with the other such operators. 


Piemulrsctestave.; witch 1s,-discissed in Chapter 2, 
is associated with symmetry operators. Because the symmetry 
operators commute with the transition operator X, as 
Welleas With H, every transition must connect states ain 
the same symmetry submatrix. Thus each symmetry submatrix 


Grmncan, be treated as, a separate calculation. 


The remaining two stages of factoring, which are 
associated with two kinds of total spin operator, are 
discussed in Chapter 3, With a suitable choice of basis, 
they can be applied_as readily to each symmetry submatrix 
as to the complete Hamiltonian matrix of a molecule without 


Symmetry. 


The rae result,of.the factoring.is, that 
ausaneleylarge,agdifficult, calculation,,involvingwthe 
complete Hamiltonian matrix, is replaced by a series 
of smaller, easier calculations, involving the submatrix 
blocks along the main diagonal. In the example already 
used, the perfluoro-2-butenes, no block larger than 4x4 
need. be; diagonalized, after all three stages. of. factoring 
have been applied. The total number of Hamiltonian matrix 
elements to be evaluated is 32896 without factoring, and 


168 with full use of all three stages of factoring. 
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CHAPTER, « Zr. SYMMETRY FACTORING 


This Chapter presents a new method of calculating 
the NMR spectra of all molecules with twofold symmetry. 
The spin Hamiltonian used is valid for all combinations 
of magnetic equivalence and twofold symmetry, including 
symmetrically equivalent pairs of magnetically equivalent 
groups. Expressions for all Hamiltonian matrix elements 
and transition intensities are given in a form suitable 


for computer calculations. 


Dak. The NMR Symmetry Group 


The symmetry properties of the NMR Hamiltonian 
were first used to simplify calculation of the NMR spectra 
Ot small, rigid molecules, such as te al co oe ee ee 
The symmetry of these molecules can be adequately represented 
by a conventional point group (C. for 1,1-difluoroethylene), 
but there is usually no suitable point group for nonrigid 
molecules such as those with free rotation about a skeletal 


single bond. 


Longuet -Higgins” has developed a more general 
approach to symmetry to describe the infrared spectra 
of nonrigid molecules, and Woodman® has shown that it 


can readily be adapted to NMR. Woodman's NMR Symmetry 


es 7 > ee Vn) We — . =~ 
oo = ®) 3 a Ray sh 7 


>a 0 F : o / a 
— 7 - t : oO a aa’ | : 
: 7 ” si a4 - 

: 2 oh 7 p 7 7 

: : ‘ 
- ——— : ra ; rte: ; 
bd “= sy a ‘ ° oy 
t r é 1 ty 


if | vn 
gist oust ns ‘to. bet pit: WSR \G einege Buy 0,83 
3 ers4emee Bioko). AoE anti i 20 LaMe pip 
saisdenhdno2 rik rol bLIGY es eau 
gated: t r ¥ rsetanes biotiows ant oon 
lmolayiupe (J laoitonash 79 eriee notte 


2themeda ziToem naruosl LNB Uf “ ok ator ee 


7 


etdetive mWrot ss. nie OWED $1.8 “aia wai 
a ayiS inetvande 4 


dew des L taal fii sits tq 208 Rates eal 


‘ ; px 
oyicoqa’ PH odd) To rev rsigalsa Yrebqoie es: Beals 


L zz : of § : ‘ — ~ 
,4ae sition out a* t,t 25 cepa cea 


t 


“hermozonqer 4h: sinupabr. ud ae arte ety bodes 


x (nina ees aout Hit > 20%: gD" quite: Fito. ‘ 
bigisuGe’ xo? quote LOY ahaa tire, Gm nn 
- feteine & ivods mp hatay- 993 Arter. a dla! 
7 ‘ " | a “t j : r: : otal 7 
hh ale sx c beqataved xi 
re 
pws 
AR 
‘e - 
_ : i 
- i _ . 


10 


Group (NMRSG) contains all permutations, of isotopically 
identical magnetic nuclei, that do not alter chemical 
shifts or coupling constants. The symmetry (or lack of 
symmetry) of non-magnetic nuclei is relevant only to 

the extent that it affects the environment of the magnetic 


nuclei. 


The correct NMRSG for a nonrigid molecule cannot 
always be predicted from the structural formula. Many 
chemical shifts and coupling constants vary with the 
conformation of the molecule. If the molecule changes 
conformation rapidly,’ the chemical shifts and coupling 
constants for the observed spectrum will be population- 
weighted averages over the values for the individual 
conformations. As a result, the NMRSG will often contain 
Syumerry clements that arecnot present ini any single 
molecular conformation. When the populations of the 
different conformations are not known, the NMRSG can 
be determined only by comparing the observed spectrum 


with trial spectra calculated for various possible symmetries. 


t A nonrigid molecule's conformation is changing rapidly if 


its spectrum shows no effects that depend on the average 
iMtetine of the conformations. In the other extreme, a 
conformation whose lifetime is much longer than the lifetime 
Siva nuclear spin state 1s in effect a rigid molecule with 
its own spectrum. Conformation lifetimes between these 
extremes are beyond the scope of this thesis; see, for 


example, pss, Chanten 0; 
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The symmetry operations of the NMRSG can be separated 
into two subgroups, either or both of which may include 
only the identity. The first of these subgroups includes 
only permutations within fully equivalent groups! of 
nuclei. The nuclei in a fully equivalent group must 
have: 

iecdual chentcalysii tts, 

ii) equal coupling constants to any nucleus outside 

Cheroroup, 
iii) equal coupling constants Within ytheveroup. 

When these three conditions are satisfied, the Hamiltonian 
Pemiiivariant to any permutation of the fully equivalent 
mucte.. che subgroup of the NMRSG that describes. the 
local symmetry of fully equivalent groups of nuclei can 
then be written as the direct product of symmetric’ 


permutation groups, ie 


Mie wturstatwo condi tions.1 oT stil lecouiva lence 
are the necessary and sufficient conditions for the 
Gonservation Oo the square of the total spin’ Of une 


group of nuclei. Every fully equivalent group is thus, 


8,9 


by definition, a magnetically equivalent group. ile 


NR ee ed 
t The word ‘group' is used throughout this thesis both for 
collections of nuclei and for groups (in the mathematical 
sense) of symmetry operations. It should be obvious from 


the context which meaning applies. 


: See §2.2 and Appendix A. 
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can be shown (see Corio?” p. 1/75) that coupling constants 
within a magnetically equivalent group have no effect 

on transition frequencies or intensities. “It will there- 
fore be assumed throughout this thesis that all coupling 
constants within a magnetically equivalent group are 
effectively zero, and no distinction will be made between 
a magnetically equivalent group and a fully equivalent 
Proupsot tnuclei. ofThe -NMRSG ott athis ythesis - isirthe 


same as Woodman's Effective NMR Symmetry Group.’ 


The second ('frame'’) subgroup of the NMRSG 
permutes “entire magnetically equivalent ‘groups of nuclei. 
Magnetically equivalent groups interchanged by frame 
symmetry operations form a symmetrically equivalent set. 
symmetrically equivalent sets of single nuclei are 
included in this definition, since a single nucleus can 
be considered a magnetically equivalent group of one. 

For example, in the perfluoro-2-butenes the single element 
of frame symmetry interchanges the two CF. groups and 
simultaneously interchanges the two olefinic fluorines. 
The frame subgroup is normally isomorphic with one of 


the geometrical point groups (C., C or C; in.the example). 


is 
This makes it convenient to borrow the notation of an 
appropriate geometrical point group, but it should not 
be assumed that any likely conformation has the geometry 


appropriate to the point group used. All that can be 


said is that the average values of the chemical shifts 


pd 


mN 


h— 
toon 
7 
** 
~ 
—t 
’ 
= 
+ 
- 
» 


-gteild 
: a 


hae gone ae 


preg Li eis if gn F2ehgenesonst sii . 


onitodes Dik 3605. .226900 2254 


ae he oe sA oT Liw roy - Sits Th oe et pris,  OROE rt ~~ ot 
. Lied, 3 ow DESERT rere 
7 t i rp 
° t 2 at ere te ea et” pet “A, ao qe 


x YT tier) ers AM mel Behe, *aenboot, 264 
° ; ‘ ' : | : > ef ah ae 


a 


wd 


TORR te «uy iyi ¢ ' ompye’) Sapoee ee 
“2 : ; " 

. ~ *. = eats oka yar " 

So agHo’ tcoobey ty v PIAS Se BR ae rid yx ant" 


¥¢ 


.3oe ‘no PRV. Li eae » at i < LTS 099 es 


» 


ep pal ; e r > p 2082 08 PowEn pa “lh ; 
7 a | e* ¢ 


Li 
* 2 v : ‘ F * “ a env" at 
bee engl, ¥] eit ¢ 7 | ; Te > SEs =f Dede E> i Ge ; ‘2 nt 


" rh — -9 a . a ‘ 
: bie -dqvox: 77 owl oft eee CTP aria 
) a , \ 4 ae 
2eatsouta areairtato os sir segs 1odind ¥le2 ones E 
7 ; — B/ = ’ Ps 7 at ; ' me: bom ‘ 
to Sno azviw siilytotoet yl laagfon £f rey: ut 
‘ (elquexe, os3 ni at, ig? - ee } equery tee Vong te: : a 
i - ' ; - fc: 
gc to nortaton oa wortod os 164 ino WROD. 
re | ; ; s aa 
: ton biwore ri tad ,Guorg talon Leo rts 


a 

— - erromireg au end HOFIAMTOT HOS isan, “ 

© coe 2 a 
: a2 feds PIA Been. quer” — 


x ial * oe ye porn ‘ 


2 


and coupling constants make the Hamiltonian invariant 


tO, particular permutations of nuclei. 


Although the NMRSG of a molecule can usually best 
be described in terms of permutations, other approaches 
Mayeteac to simpler calculations. “For example, it 15 
Derreculy Valid to classify the spin’ functions of a 
magnetically equivalent group according to the irreducible 
representations of the appropriate symmetric permutation 
group oe Since the character’tables of all symmetric 
permutation groups useful in NMR are readily available,!9 
Pie pincioles there 1s nosdifticulty in. constructing a 
Gasis With the proper symmetry. ~However, in order to 
calculate Hamiltonian matrix elements and transition 
Teer atures. elt iS necessary to-expand each basis function 
dimeceris Of the’ origvinal basis-(1.2)" For’the larger 
Magnetically equivalent groups this would involve a vast 
amount of unnecessary calculation. The composite particle 
method, described in §2.2, provides identical results 


Without requiring any explicit expansion of the symmetrized 


basis. 
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2.2 Magnetic Equivalence 


In the composite particle approach to magnetic® 


equivalence, a total spin operator F is defined for each 


magnetically equivalent group: 


Se Py ten Goer) 
ah) alg Rewe 
Any single nuclei are considered to be magnetically equivalent 
groups of one, with i = Te PhesspinyHaniiconian, (is 1) 


is rewritten in terms of these total Spins Operacors: 


H = Lv 2, (2) + cVitesh(t) Els), (2.2) 
The conditions for magnetic equivalence are the necessary 
and sufficient conditions for F* to commute with H and 
Mipvechentransition operator X of Eq. .(1.3). -As a result, 

it can be shown that the terms involving coupling constants 
within magnetically equivalent groups, which have been 
Omitted from the Hamiltonian (2.2), can have no effect 


on the high-resolution NMR spectrum. 


Comparing*{1V1) and?(2.2)5*iteis apparent that 


a magnetically equivalent group, when in a state with 


T This section presents only a brief summary of the useful 


results of the composite particle method. For a more 


detailed account with proofs, see 85.4 of Corio” 
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span F, can be treated as a composite particle of spin 
fe] Fs, A magnetically equivalent group of n nuclei of 
spin I takes on all non-negative F values in the sequence 


Peele leat l= 2, 


To take advantage of the commutation properties 
of the F? operators, the basis for the system should 
Demconstructed from their eigenfunctions. A simple product 


focimeecioilear to .( 1.2) cis isuitable: 


E |Fmy?. CS 
When these basis functions are correlated with the 
Dnoyimetrized basis functions (1.2), it, is apparent that 

for each group all F values except the highest are degenerate, 
with degeneracy Sr. An index s = 1,2,..-58p can be 
introduced; the basis vectors satisfy the orthogonality 


relationships 
<Fjomj35,/F;,m53s5> 2s. acre eye Caan 


The indices i and j are used here for a magnetically 


equivalent group in two basis functions. 


Since F? for each magnetically equivalent group 
commutes with the Hamiltonian (2.2), the Hamiltonian 
matrix can be reduced to block-diagonal form with each 


block (submatrix) labelled by a unique set of F and s values. 
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2 : dey os 
Because F° also commutes with the transition operator X, 


every transition must connect energy levels in the same 


block. 


Ine > Hamiltonian (2.2). 1S written in terms, of i and Lee 
operators only; thus blocks whose basis functions have 
identical F and m values, but different s values, must yield 
identical Hamiltonian matrix elements. Similarly, the 
Meotewerois. if these blocks: must; hbevadenticalvs» The spectrum 
corresponding to each unique set of F values need be 
eeleulated"only* once); and its™antensitiessmultipliediby the 
statistical weight Igy (the index i running over magnetically 
equivalent ay he : which can be calculated without 
reference to the original |I,m> basis (see Corio” Dare LZ Uv). 
whese.partial spectra forgall sets,of.F values.are then 


superimposed to give the complete spectrum of the spin system. 


The same total spin properties of a magnetically 
equivalent group that provide useful factoring of the 
Hamiltonian matrix also lead to the new form (2.2) for 
the Hamiltonian. Matrix elements for an entire magnetically 
equivalent group can be found as easily as those of a 
single nucleus. In contrast, the permutation symmetry 
of the same magnetically equivalent group does not lead 
to a usefully simplified Hamiltonian. Instead, the original 
Hamiltonian (1.1) must be used to evaluate the complicated 


expressions for matrix elements in the symmetrized basis. 
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As the two approaches yield identical factoring,’ 
the composite particle method will be used throughout 
the remainder of this thesis. A particular set of F 
values for all magnetically equivalent groups will be 
treated as defining a unique block (submatrix) of the 
Hamiltonian, with the implicit assumption that an 
appropriate statistical weight must be applied to all 


GaeitSe transition intensities. 


C. W. Haigh has proposed (private communication) 
an excellent notation for symmetric spin systems, which 


Pseused throughout this Thesis. 


As usual,! letterssclose;, together: in: the; alphabet 
(A,B,...) are used for magnetically equivalent groups 
with similar chemical shifts, and other letters far from 
them in the alphabet (W,X,...) are used for groups widely 
Paperaced ti chemical shitt fromthe A,B... groups. 


The number of nuclei in each magnetically equivalent group 


For magnetically equivalent groups of nuclei of spin 1/2, 
there is a one-to-one correspondence between total spin and 
Prreduciple representation in ee But for groups of 
Preece or more nuclei of spin: 1 >7%, blocks, of basis functions 
degenerate with respect to F may belong to different 
irreducible representations of ae This yields no extra 
mectoring, and iS in any event:unimportant; such groups 
are seldom significantly coupled to the other magnetic 
nuclei in the molecule. See Diehl, Harris, and Tones 


footie ther references. 
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is indicated by a subscript. Thus the spin % nuclei 


in monochloroethane form an A,X spin system. 


An n-fold element of frame symmetry is indicated 
LIS rere | 


of each symmetrically equivalent set. Thus 1,3,5-trifluoro- 


ye enclosing the notation for a typical member 
benzene’ 1s an [AX] Spin system, and c1enem Oc wEnewpermLuoro. 
e@-buULenes 1S an [AX.], Spin system. Groups invariant 

to the symmetry operation are placed outside the brackets. 


Thus p-fluorotoluene could be described as Az[KL],X. 


The extension of the notation to molecules with 
more than one element of frame symmetry involves the 
use of more than one set of brackets. It will not be 
described here, because such molecules are outside the 


seunewot this Thesis. 
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ze OM Twofold Frame Symmetry 


“ghee Weal ig The Spin Hamiltonian 


The composite particle method has greatly simplified 
calculations involving magnetically equivalent prouse but 
only conventional group theory has previously been used to 
deal with frame symmetry. It is shown below that for mole- 
cules with twofold frame symmetry it is useful to apply 
mua ceneral princrples of the composite particle method: 
mies Das ls. 1S Boner rieeed from eigenfunctions of group 
total spin operators, and the Hamiltonian is written in 
terms of spin operators that embody all the symmetry of 


Eieo spin system. 


When there is twofold frame symmetry, the Hamiltonian 


(2.2) can be rearranged as follows 
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In (2.5), magnetically equivalent groups which are transformed 
into themselves by the frame symmetry (invariant groups) 
Ptemucvebicd i,),.... All: other cours must appear in 

Dore labelled p,q,..., Such that the frame symmetry 

element interchanges the primed and unprimed members 

Seer ciepdir. iwe two distinct coupling constants (usually 
Gas and trans) connecting pairs p and q are written J 


Pq 
and J... The coupling constants connecting the two members 


qP 
Creaspair »to any invariant group must be equal or the 
wumectyeis destroyed. As a result it 1s pointless» to 
becewirome Symmetry unless there are at least two pairs 
momaccsit1on to-any invariant groups. A _single pair, 
as in 1,3-dichloro-benzene, meets the conditions for 


Magnetic equivalence and is better treated by the usual 


composite particle method. 


The actual nature of the twofold symmetry operation 
is irrelevant as long as the Hamiltonian can be written 
in the symmetrical form (2.5). In practice the NMR symmetry 
is usually the result of a plane or twofold axis of symmetry 
in the molecule, taking the conformation as an average 
over an NMR time scale. For convenience, the notation 
of the group C. will be used below, but this should not 


be considered a restriction to molecules with that symmetry. 


7A cCOLalL spin operator A is defined for each 


symmetrically equivalent pair of groups 


A(p) = F(p) +F(p’), (2.6) 
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and a new basis is constructed which diagonalizes A? 
lL 
and A, for each pair.’ In terms of the basis (2.3), 


the new basis for each pair may be written 


[A,M> = } (F,F*,«,M-«|F,F*,A,M) |F,«>|]F°,M-«>. (2.7a) 
K 
The Clebsch-Gordan coefficients’ (and the |A,M>) exist 


only for values of A satisfying the triangle inequalities 
|F-F*| < A <$ F+E7 (276) 


The basis for the entire spin system then consists of 
a product of magnetic equivalence functions for all 
invariant groups and pair total spin functions for all 


symmetrically equivalent (pairs ‘of ‘groups 
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Bits a 


t It should be emphasized that, by definition, none of the 


A? (p) commute with the Hamiltonian. If A* for any pair of 
Broups does commute withsi, allethe nuclei;sin the pair of 


groups form a single magnetically equivalent group. 


Y The notation of Condon and Shortley!? is used for the 


Clebsch-Gordan coefficients (F,F*,«,M-«|F,F~*,A,M). 
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One can now rewrite the Hamiltonian (2.5) using the symmetric 
spin operators A(p) and the complementary antisymmetric 


operators 

B(p)) =» Fp) -- F@*), (2.9) 
obtaining 

Ht, = WyiE@) + aged EG) 
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+ J Jo, FG) -A@) 
ap 


eee) } Typ fh) EO) (2.10) 


+% JY (J. 43. JA(p)+ACq) 
p<q PQ qp 
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The transition operator is now written 
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ee. Matrix Elements in the Symmetrized Basis 


Many of the matrix elements appearing in (2.10) 


mre casy tO é€valuate. FE. and A, can contribute only 


to diagonal matrix elements, since the basis is constructed 


from their eigenfunctions. FF, which commutes with 


both A? and A,, also contributes only diagonal matrix 


elements, whe first formulas 12%2 of Condonfand Shortley +? 


becomes, in the present notation, 
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<A,M|F+F°|A,M> = %[A(A+1)-F(EF+1)-F*(F*+1)]. (2eaa) 


Pnesremaimming -terms-an (2.510) can? contribute’ to 
either diagonal or off-diagonal matrix elements. Their 


POEsproOduuct, tOrm permits the separate evaluation of the 


two halves. of: each product.» Each half\issa matrix) element 


such as <A,M|A|A‘,M*>, whose vector character is shown 
in the following equations by the use of the Cartesian 


, + 3 > 
Diiieveoctors 1, 3),k. 


The matrix*elements of F for the invariant groups 
are well known. Those of A and B for the pairs canjbe 
found from the equations for matrix elements of Ff and F- 
in an |A,M> basis, given in pages 56-73 of Condon and 
Shortley.!2 The matrix elements of A in the |A,M> basis 
are, of course, similar to those of F in the |F,m> basis. 
Since A commutes with A’, AA=0 for all its. nonzero matrix 


elements: 
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<A,M|A|A,M> = MK 


‘Coats 
<A,MJA|A,Mf1> = [ (AM) (AM+1)]® Z#23. 
ssl 


The commutation rules for B with TeSpeECr, to A follow the 


common pattern 


[A, B,J = 0 [A, sBy} = 4B Md) aera 
eae 0S ate [See Sle, Seeamue re, eam 
ee aie Based 1 fe [A By] i a 

Pema result, 
<A,M|B| ALAA ,MtAM> = 0 unless AA=0,1;AM=0,1. (2.14) 


Nonzero matrix elements of B allowed by (2.14) are 


1 
<A,M|B|A-1,M> = §(F,F*,A)[A?-M?]* k 


<A,M|B|A-1,M#1> {(F,F7,A) [ (AEM) (AFM-1)] ® 2243 
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1 
<A ,M|B|A,M+1> g (F,F~,A) [ (AFM) (AM+1)]* T#2j 
~ 2 


<A,M|B|A,M> = g(F,F*,A)M i 
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where 


[A2-(F-F7)?] [(B+F7+1)?-a2] |7 
A? (4A2-1) 


(Ge A) 


C2210) 
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A(A+1) 


g(F,F°,A) 


Note that when F=F”, g(F,F°,A) vanishes, giving the stronger 
condition AA=+1 for nonzero matrix elements of B. The 


expression for {(F,F°,A) then reduces to the simpler form 


§(F,F°,A) = | (2F+1)?-a? |” (2.17) 
4A? -] 


As an aid to hand calculations, the squares of the factors 
ine (2.16) Gre given in Table 2Z-L for all combinations 
Cue, F ,4end A that occur for pairs of@groups with spin 
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Squares of Factors! ACB Re Ay) -and. 9 CE shes) 


of <A,M|B|A~,M*> 


i See equations (2.15) and (2.16). Primes on A and M 


are used only to distinguish between bra and ket values. 
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Sit Wiehe wr Symmetry Properties of the [| A, M> Basis 


When the molecular twofold Symmetry operator C, 
is applied to the basis (2.8), the basis function C274 


for each Symmetrically equivalent pair becomes 


C, |A,M> = }(F,F*,«,M-«|F,F*,A,M) |F*,M-«>|F,«> (Zars) 
K 


The resulting basis function clearly belongs to a submatrix 
of the Hamiltonian related. to the original submatrix by the 
interchange.of all primed .and »unprimed ;groups.,:,There are 
two possible cases: when for every pair F=F’, the new and 
old basis functions are the same; otherwise, C, relates 


hasassfunctaons «in.diféferent submatrices. 


Thus when F#F°~ for at least one pair, the product 
functions (2.8) are not eigenfunctions of C,, and the 
magnetic equivalence submatrix cannot be further factored 
on symmetry. However, from the symmetry of (2.5) it is 
obvious that the two submatrices related by C, are 
degenerate. Energy levels and transition intensities need 


be found for only one member of each such pair of 


submatrices, giving it double the usual statistical weight. 


When for all pairs F=F”, the redundant prime in (2.18) 
can be dropped. Substituting m=M-k and reversing the order 


of summation, one obtains from (2.18) 


C, |A,M> = )(F,F,M-m,m|F,F,A,M)|F,m>|F,M-m> (2.19) 
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Interchanging the two F values, along with the corresponding 
m values, at most changes the sign of a Clebsch-Gordan 


coefficient.' Thus (2.19) becomes 


Cy |A,M> = J¢-1)7F Ace, F,m,M-m|F,F,A,M) [F,m>|F,M-m> 
m 
(2.20) 


GENE 


(1s) 
For each pair, |A,M> belongs to one of the irreducible 
representations of C.. The state of highest spin belongs 
to the A representation, and those of lower spin are 
alternately B and A. The complete product function (2.8) 
will have B symmetry if an odd number of the individual pair 
functions have B symmetry; otherwise it will have A 
Symmetry. The invariant groups, which always have A symmetry, 


cannot affect the symmetry of the product function. 


Every submatrix in which all F=F” may now be 
separated into two submatrices containing the symmetric (A) 
and antisymmetric (B) basis functions respectively. Because 
C, commutes with both the Hamiltonian and the transition 
operator, these submatrices are completely independent of 
one another, just as they are independent 


of all magnetic equivalence submatrices with different F 


values. 
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20504. Transitions 


Transition intensities are calculated from the 
matrix elements of the transition operator X (see Eq. 2.11). 
Transitions involving invariant groups have intensities 


proportional to 


| Spe U Ea rem- leet? [oe a) 


and those involving symmetrically equivalent pairs have 


Belensutaes proportional to 
| <A,M|A, |A,M-1> ie (22159 


All F* commute with X, giving the selection rule 
AF=0. Similarly, for those magnetic equivalence submatrices 
inewarcn all F=F~ , there is: the additional selection) rule 
that the symmetry must not change. Both of these selection 
Tiveseare taken care of implicitly by treating each 
submatrix ds “an,independent calculation. It must be 
remembered that all transition intensities for any submatrix 
are multiplied by the appropriate statistical weight to 


account for any degeneracy. 


Since; all A* commute with X, there 1s the additional 
selection rule within each submatrix that AA=0 for every 


symmetrically equivalent pair. 
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The=torm of the tsymnetrrzed tamiltonian, Eq. (2.10), 
was chosen to show its derivation clearly, and as a result 
it appears unnecessarily complicated. It can be simplified 
Pyeredeftining -some of its terms. The definition of the 
Poca sin: Operator A can be extended to include F for 
Picwonvariant groups, just as id includes i forall sing te 
Reaolot es Ouls Over 1.,j)... now include» both invariant 
Proups, and symmetrically equivalent pairs; sums restric- 
weuetor ric pairs now Use the indices 2,m,.... in each 
case a symmetrically equivalent pair is considered to 
be a single entity, represented by a single index; primed 
indices such as £° ,m° are used only where the two halves 
of the pair must be distinguished. The two coupling 


constants Jog and J are combined to give the symmetrized 


qp 
couplings 
eet 
Jin = fs tpaignelap? 
C22.) 
=- 1 ~ 
Jme = #9 pq ~ Jap? 


With these changes, the symmetrized Hamiltonian (25.0) eoan 


be replaced by 
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iieebasis (2.8) reduces to 


| A. ,M,>, (2:24) 
1 
and (2.2la) is no longer needed. It must be remembered 


that the two F values for each symmetrically equivalent 
Davo tenessential to thesdefinition of the basis, and 


appear in the equations for many of the matrix elements. 


Equation (2.23) shows more clearly the distinction 
between a symmetrically equivalent pair and an invariant 
Pours ALi Matrix elements of A phy eles Bo a) aly ete seems a Key g ea oe: We 
to those that would be. found for an invariant group 
Peres mrie ofall lathe nucle: in the-pair.. The coupling 
constant J om that appears in A matrix elements connecting 
twO Pairs is just the average of the near and far couplings. 
But since the difference coupling Jnl is not zero, matrix 
eLements. of B can mix basis functions where a pair has 
different A values. This in turn allows the coupling 


Peeewurninathe pair to aftece the appearance of the 
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spectrum. 
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Frame symmetry does not factor all of the magnetic 
equivalence submatrices when there are symmetrically 
equivalent pairs of magnetically equivalent groups. However 
the largest submatrix, the one in which every group has 
its maximum F value, can always be separated into symmetric 
and antisymmetric parts. In these and all other symmetry- 


7 


Peclored submatrices, F =" FS for every*pair. AS a result, 
B contributes to off-diagonal matrix elements (AA = #1) 
only, and the simpler equation (2.17) can always be used 
Pwececulate.them.. If «v.25 )easgtused for the submatrices 
which are not factored by symmetry, the full equations 
ize. pnd.({2.10).must be.used forsat .deastjsome B matrix 
elements. No factoring is lost here if the unsymmetrized 
Hamiltonian (2.5) is used, but the advantage of simpler 


equations for the Hamiltonian matrix elements may be 


outweighed by the need to treat each pair as two groups. 


The usefulness of the symmetrized Hamiltonian 
and basis derived in this Chapter, is significantly 
increased by its suitability for computer calculations. 
Two characteristics of the method are particularly 
important in this respect. First, the simple product 
form of the basis (2.24) is valid for molecules of 
arbitrary size. Second, because the two halves of each 
dot-product Hamiltonian matrix element an (Qiao) ‘are 
evaluated separately, the calculation of each such matrix 


element depends only on the characteristics of the two 
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groups (or pairs) involved, regardless of the size of 


the molecule to which they belong. 


PeCoMputer program based on (2.23) and (2.24) 
can, therefore, calculate the spectrum of any molecule 
with twofold symmetry, subject only to limitations imposed 
Dreeouputer size. Because (2.23) and, (2.24) include 
(ene eatioe(2..9), the Same program can equally readily 


calculate the spectra of molecules that lack frame symmetry. 


ttwis possibie that) C symmetry could be treated 


ZV 
by methods analogous to those of §2.3, but the many possible 
combinations of local symmetry functions in the molecular 
product functions would lead to a very complicated set 

of equations for Hamiltonian matrix elements. For most 
Otiererrame symmetries, not all molecular product functions 
iitaistormeas one of the irreducible. representations of 
ticemolvecilar symmetry group. -Thesproduct functions 

hOmeany particular spin system, such as [AX] z, can be 
resymmetrized to yield correct symmetry functions for 

the whole molecule ,” Ditet eect ke i sthatache process 
could be reduced to an algorithm suitable for a general 
computer program. Special programs could be written, 


but few of these spin systems are common enough to warrant 
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CHAPTER 3. FULL FACTORING 


1 The Hamiltonian Matrix 


The NMR Symmetry Group provides only the first 
stage of factoring the spin Hamiltonian. To permit easy 
Peserence, to each Stage, a simple subscript ieee 
toevecdy in) the remainder of this Thesis... The complete 
Spin Hamiltonian matrix, prior to any factoring, is referred 
to as the Hamiltonian matrix, . Symmetry factoring yields 
matrices., Woiechsane factored by=the z component of fotal 
spin, M, into matrices,. When. there. are .large chemical 
shift separations, each matrix. 
Poeeeparticular form of factoring does not apply, some 
Giactiesceterms become synonymous. For example, if there 


Meeioesyimetry factoring, the single matrix, includes 


the entire matrix). 


Pte seconvenicentstosextendsthis#notetion: to include 
other quantities associated with the Hamiltonian matrix. 
Thus basis functions, eigenvalues, eigenvectors, and 
transitions are also referred to as. belonging to particular 


matrices,, matrices,, and matrices). 


Only one kind of factoring always applies: every 
matrix, ‘oLactoredsinta matrices., each labelled by a 
different eigenvalue M of the z component of molecular 


Gotdlesuis. » nere. are 2M yt matrices, in each matrix,, 


Teese ho redai nto matrices,. 
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with M = M M -1, 


+.49inM ; 
max max 


to the sum of the F values of all magnetically equivalent 
groups in the matrix,. (When there is twofold frame 
symmetry, Me in an antisymmetric matrix, is one less 


than M insthe symmetric matrix, with the same set 


ax 
of F values.) 

ane Size of Each matrix, is determined by the 
number of possible basis functions with that M value. 
When all groups have F = %, with each basis function 


one of the 2” possible combinations of a and 8, these 


Sizes) 20 a5 the appropriate binomial coefficients! “In 
Penerat it Can’ only be said thatthe matrices, with 
M = +M are smallest, and those with M nearest zero 


max 
Brestarpest,, Within each matrix,. 
whesthird stage, X factoring ,.derives its name 
from the usual convention! Or destribing nuclear spin 
systems. Nuclei forming a set with similar chemical 
Siareso ate, labelléd A,B,...; those in a second” sét are 
Dae ted. AG) ,e. 2 LL all of the coupling constants Jay 
connecting them to the first Set are much smaller than 
the corresponding chemical shift separations Any: The 
concept can be generalized to: allow any number of sets 
of nuclei, each of which meets these conditions with 


Tespect. to every other set. “EaGh such’set of nuclei” in 


a molecule is referred to here as a.species. 


The basic assumption of X factoring is that AGS); 
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the z component of total spin for each species, commutes 
with the Hamiltonian. Each corresponding eigenvalue, 

Mo; is then a good quantum number. This approximation 

is valid only when all terms containing the first and 
higher powers of Tax/ Any can + legitimately be dropped 

in a perturbation calculation (see Corio” Chapter 7a. 
PiecocevariOus Specieseareniy fact different isotropic 
Species'’such as '*F and *'P, the chemical shift separation 
is many MHz at normal field strengths, compared with 
coupling constants of at most a few KHz. The X approximation 
is then essentially exact. Within proton spectra, coupling 
constants are usually small enough to allow X factoring 
when Any is only a few hundred Hz; for other nuclei, 
particularly ‘°F, the wide range of chemical shifts makes 


X factoring even more common. 


When all the M, values are good quantum numbers, 


S 


Bae be eeunccions that differ in. any Mg value cannot mix. 


Week atactoxred Matrix, reduces) to a, series of matrices,, 


3 
each of which is characterized by a unique set of Mg 
values whose sum is M for that matrix,. Fig. 3-1 shows 
the resulting factoring of a matrix, from an [AX.], spin 
system (the antisymmetric matrix, WiC we, S504 BA % , 

Fy = Fy. = 3/2). The two My 
and the M value for each matrix, are given beside the 


values foxy each matrix, 


matrices to which they belong. The matrix, boundaries 
show the amount of factoring that would be available 


if the two chemical shifts were too close together to 
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FIG. 3-1. FULL FACTORING. 
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permit X factoring (ire. for<an [AB], Spin system). 


The matrix sizes in Fig. 3-1 illustrate several 


Pypecat points aboute®xX factoring. A large matrix, is 
generally factored into more matrices, than is a small 
matrix,, but this does not prevent an increase in matrix, 
size as the matrix, Size increases. Because some of 

the matrices, are GuLte SMaibs the largest, one_in each 
matrix, may have a dimension only slightly smaller than 
the matrix,. There is however a considerable saving 

in calculation. In the example, the three matrices, 
in the matrix, with M = 0 have sizes 1x1, 4x4, 1x1; there 
are, therefore, no more than 18 nonzero matrix elements 


compared with 36 for that matrix, without X factoring. 


3 
iii Ddsis Lunctlons. an a matrix, can be paired, 
associating each one with its complement in which the 
oe eos. all the group M’values are reversed.” “since 
these basis functions have complementary M and Mg values, 


the pattern of matrix, and matrix, Sizes.is Symmetrical 


& 
about M = 0. In hand calculations, once the Hamiltonian 
matrix elements for each matrix, have been found, those 
for the complementary matrix, can be obtained much more 
quickly. Every term will be identical in magnitude to 


the corresponding term in the complementary matrix, , 


but some terms will have the opposite sign. 


Theemost extreme % Lactoring. OCccurs,-1or a molécule 


in which no two magnetically equivalent groups have closely 
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spaced chemical shifts. Since each magnetically equivalent 
group is then a member of a different species, every 

group M is also an Mg+ No two basis functions share 

the same set of F values and the same set of Mo values, 

so the Hamiltonian matrix can have no off-diagonal elements. 
Thetrésultris:a simple spectrum of:the type given in 

Shl2 o¢ pst, Much of the popularity of NMR as an analytical 
Pooieicsuits trom, the fact that the marcnitude. of, every 
COuUpiadeeconstant canbe taken directly from the separation 


of experimental lines when a spectrum meets these first 


ay 2 
order conditions. 


The spectrum of a molecule with frame symmetry 
is less easily analyzed, even when there is only one 
dastinct)chemicalashiftnper spéciestheThextwo or more 
magnetically equivalent groups in a symmetrically 
equivalent set have identical chemical shifts, and are 
therefore of the same species. Thus even the simplest 
of all spin systems with frame symmetry, [AX], ; never 
yields a true first order spectrum. 


One additional type of factoring, which has been 


applied by Harris!4 


to some [AX], spin systems, should 

be mentioned briefly. In strongly X factored spin systems, 
it may be possible to separate all the basis functions 

of a matrix, into groups, such that every off-diagonal 


matrix element between basis functions in different groups 


is some multiple of one coupling constant. If that coupling 
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constant is zero, each of the groups of basis functions 
becomes a separate matYrix,. In hand calculations this 
extra factoring can yield additional easily solved 1x1 
and 2x2 matrices. In computer calculations it is best 
ignored for two reasons: the matrices, that can be further 
factored are usually quite small, and there is no simple 
rule for sorting the basis functions into matrices,;. 
Zero coupling constants produce additional factoring 
Only @iE they appear in off-diagonal’ matrix elements’. 
The basis that allows this may be different from the one 
that would normally be used when all coupling constants 
are nonzero. In the twofold frame symmetry basis of 
Ghapter.2, both Jane and Jyy- appear only in diagonal 
matrix elements, and cannot affect the degree of factoring. 
There are 3x3 and 4x4 matrices, to diagonalize. in an 


[AX sy stem. .and--Larger,ones fom Langer [AX,], systems. 


3]2 
However if a simple magnetic equivalence basis is used 
for any [AX], system with Jyy- = 0, there is no matrix, 
larger than 2x2. When Jyy- is nonzero but small, as 

: 14 . 
in some molecules discussed in 85.3, Harris's algebraic 
solution for Jyy- =.0..can, still.be. applied, to, obtain 


a first approximation to the correct parameters. 
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It was pointed out in Chapter 2 that the symmetry 
factoring of the Hamiltonian matrix is exactly matched 
by transition selection rules:, each matrix, 1S.a separate 
entity whose spectrum can be calculated without reference 
Poet hesocney, Matrices ,.«, Within each matrix,, the factoring 


2 


by M and Mo into matrices, and matrices, is also reflected 


3 
in selection rules. These differ from the symmetry 
Seeew tO rules, chiefly insthatatransitions occur, not 


within a single matrix. or matrix,, but rather between 


S 


certain matrices, in adjacent matrices ,. It is. shown 
below that these selection rules can be derived quite 


Sie ivetrom the nature of the-factoring described in 83.1. 


The frequency of a transition is easily obtained 
from the algebraic difference in energy of the pair of 
states the transition connects. The corresponding intensity 
is obtained from the absolute square of the matrix element 
Die Weausee dec i. ZI) ed between these two states. This 
matrix element is best found indirectly. First, the 
readily calculated matrix elements of X in the original 
basis. (2.3) are obtained... Each; desired matrix element 
is then calculated from this initial representation of 
X using the unitary transformation that diagonalizes 


the Hamiltonian. 


The transition selection rules in the basis representation 
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are straightforward. There can be a nonzero matrix element 
of X between a pair of basis functions only if? 


i) AF = 0 for all magnetically equivalent groups 


il 


11) AA 0 for all symmetrically equivalent pairs 
-&id ) AM; eet lLatonrsomesones (invariant) sanqup,or 
(symmetrically equivalent) pair 
iv) es selatorealleother eroupseand, pairs, 
Since, in the basis representation, the .Hamiltonian:matrix 
consists of matrix, blocks along the main diagonal, the 
unitary transformation that diagonalizes it can mix 


basis functions only within a matrix,. Thusvenesesselection 


Mises soperate. at the matrix, beveleatter then trans: ormation, 


Basis functions that meet the first two conditions 
must be in the same matrix,. To .meet.the, remaining two 
conditions, they must be in different matrices, : 1 
DVeeetiecOrhcroupei, but snot for anysother greup,inythe 
Same species as group i, then the Mg Vinvece fore thatespecaos 
G@icereby 2lo4thesfourth condition, requires. that AM, = 0 
for every species other than the one containing group 
i. After transformation to the representation in which 
the Hamiltonian matrix is diagonal, it may no longer 
be possible to label any particular transition as belonging 
Poeeroupe 1, .but all ofythe transitions connecting, the 
Daix of matrices, must still belong,to the species whose 


Mo changes. 


Matrices, that meet these conditions (AM, = +] 
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fon One especies , AM, = 0 for all other species) are found 
in matrices, whose M values differ by +1. If X factoring 
iS net wsed; Mo for the "onlys pec tes tis “ident ical<to 
M, and an argument similar to the one above leads to 


1 


the usual selection rule AM = +1.. Transitions can then 


connect any pair of energy levels in adjacent matrices,. 
Thus the additional selection rules introduced by the 


X approximation are a natural consequence of the additional 


factoring of the Hamiltonian matrix. 


Dtewase DOInted Out Unessnlethat x factoring. 1s 
an approximate method in which all terms in Jay/Any are 
assumed to be negligibly small. As in any perturbation 
method, eigenvectors are affected more than eigenvalues 
when the neglected terms begin to become important. In 
an NMR calculation, the eigenvectors are used only to 
transform the matrix elements of X (see Eq. (2.21) ) 
to the representation in which the Hamiltonian is diagonal. 
Thus fea calculation “bs:-carried “out wsing X factoring 
tiateicenot quite justified by the chemical shifts and 
coupling constants, the intensities of the calculated 
lines may be somewhat distorted, but the frequencies 
will be very close to the true values. As a result, 
if the observed and calculated spectra are compared chiefly 
with respect to frequency, X factoring can be used to 
save time with little sacrifice in accuracy. A good example 


of the use of such borderline X factoring is given in 
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[eeeiiecberalalysis of the spectra, of the 2,5-dimethyl- 


3-hexenes. 


If the X-approximation is valid for a molecule, 
the transitions of each species are grouped together, 
far from the transitions of any other species. In 
calculations it is usually convenient to follow experimental 
practice, and give all chemical shifts for a species 
as displacements from a reference frequency in that species' 
Terpson' of the spectrum. Usually the reference for a 
calculation will be the same one used experimentally, 
DitereeisesOMmetimes easier to use one of the: chemical 


Saattes in the molecute instead. 


The chemical shift contribution to the Hamiltonian 
Wee oiecan then be written in terms ‘of the species reference 
frequencies Vo and displacements V5 of each chemical 


shift from Vo Ota Ss! 5 peGles. 


IvjA, (4) = TvgA, (8S) + TVA, (4). (3.1) 
i S 1 

Theefaretesum om the right of (3-1), which contributes 

equally to every diagonal element of a matrix,, contains 

all of the information about the actual separation between 

species. The selection rules described below allow 

transitions of species S only between matrices, where 

this Sum differs by tV3 thus each species' spectrum is 


unaffected by the vo values of the other species. Usually 


44 


be o | 
-[qtsaeth-2,§, od to at? 33% odd age =m 


‘ e | i. & ow » * La 
etpasier a tot bitey . 
| 4 : > a) ‘ ru T ° 4 - —§ p 
TonTsgos: Vege’: * ‘ 
jet ceaboaq? rari | 


ae 7 - ne igne ae i. 
f se vitwueureh ai efeotvaiae 


fatasmiTsqrs worse aa ee “Ute eel os 


+ 
>’ 
- 


’ ee ® ; = Pe a = 
97904 , por ite Lusrate "ae ove a bee ‘ ve : 
dee St OS a? : on : 
= ' 7. ory re ? y 3 |? 
Teatnaqa jalld me one tt radites a maw 2 is09 76 be 
: + ox “| yoerrege ade te = 
L rol , 5 j : Ae 
'tetasmiredro Wee ws anal one Sek Shi waeees 
4 ‘ 2a oe t ‘ | 
. of sHiess Sent eSRee eee 


foxitoedo 3 1 9 0 : 
a Ga 4 = [os 


beol uti’ sfuogion ms. ak 2 


¥ : 
hae i 1 f Lavoe 
fet Oy Loe, 

t - 
.* are ent 2 7” ie tate 3 r ; . |. ff BES [ LN 

IMTS ToT ef 9994 ; ey 

F : Die 
( Lentureno- ao8e , | a4 : 
cande aa) aa ay 
’ Hi 
— — ’ . 7 s 7) ~*h oe | —_ 
; 5 7 j 3 bs ' < = ' “a sang? of 
| f ae 


a a% se 7 , Joe — 
sesudivados dain , (1-6) 1 iigia sit 1g me a 


ziti nlne> ey eht io &-1o snehels LenQogee ED: a i sd 


7 % i, oi a” 
a2 [Rotor oat; tulos is no nie i, 
“. bi 
° i 


‘a ‘i wotte woled bedttozeb eu lut aapeanaal 
tal be | : 


neswsed qobte ru 


; z i —_: am 
a Totty atta nacweod yYiao ¢ 


oy F age Asso 2utht ‘pat 
| i ek mizzeq2 * ‘gorau "4% #9 ss igve 
7 at ing ae "a - 
4 t tama + o0q2 revo ot to 2st AV 
P. nent ee) en ees | sais in, 
S en a : >} , a 


45 


all Vo Values are set ‘to.zero,. and Va values are used 
bore cneachemical shifts, \Aliecaliculated transition 
frequencies of a species are then displacements from 


Vo rather than absolute transition frequencies. 
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be Ppa Phe olLructure OL Lig. opectrum 


pre er Bes Subspectra 


Together with symmetry factoring, X factoring 
breaks up the complete spectrum into smaller and more 
manageable groups of transitions. A particular matrix, 
yields a set of transitions that can be calculated without 
reference to. the other matrices,. When. there, is. X, factoring, 
these transitions can be separated according to species. 

Any single species' transitions can be further separated 
ijeoeseis, Cach characterized by -a.particular. set. of Mg 
Moose Owl other species..o lhe factoring thus, imposes 
aestructure on, the spectrum,, making, certain, lines. (not 


necessarily close together in frequency) members of natural 


groupings. 


Each such natural grouping of lines, characterized 
by a particular set of quantum numbers, will be referred 
to as a subspectrum. ++ Though most references to subspectra 
Wilbepemto. smaller oroupings: of lines... all of the lines 
belonging to a particular matrix, could also be called a 
subspectrum within this definition. The structure of 
the complete spectrum in terms of subspectra is independent 
of the chemical shifts and coupling constants, as long 


seni nerconditions.or Lactorning, are. preserved. 
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Because a subspectrum is defined by a set of constant 
quantum numbers, its pattern of lines may be strongly 
aifected by only a few of the chemicalsshifts and coupling 
constants. Often the functional relationship between these 
parameters and the line positions is simple enough to be 
direc#ly useful in the analysis of complicated spectra. 


Many of the methods of Subspectral Analysis?! 


are based 
on similarities between particular subspectra and the 


complete spectra of well known simple spin systems. 


Fig. 3-2 shows thegpattern of allaqwed transitions 
in the symmetric matrix, with all F = % from an [AX], 
system. This matrix,, whoseestatistL.calaweicht as four ; 
aseidentical in structure. to the symmetric matrix, of 
an [AX], Suimesystem. -Ihe expressions. fer transi tion 
Preduencies and intensities in an [AX], spin system, 
Pavengon p. 141 of pspt , are equally valid here. If 
the transitions of this subspectrum can be recognized 
gmeecithner the species A- spectrum or the species X spectrum 
of the [AX], system, all of the coupling constant combinations 


except Lane -Jyx-| can be calculated. 


The species A transitions of this same matrix, 
are shown again in fig. 3-3(a), separated into subspectra. 
The simplest subspectra correspond to the extreme values 
for My 5 this is generally true in larger systems as well, 
partly because there are fewer basis functions with extreme 
Mo values than with Mo Valueseicare7ero-.. Ln chts. case, 
the lines with My = +1 form a strong, readily recognized 
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SPECIES A TRANSITIONS 


[AX 1, SYMMETRIC MATRIX, ALL F=1/2, WEIGHT 4. 


FIG. 3-2. ALLOWED TRANSITIONS. 
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doublet. The remaining four lines, whose total intensity 
is half that of the strong doublet, form the My = 0 
subspectrum. 

Within a particular species, the matrices), as 
well as the-transitions,.can be labelled _as_belonging 
POseubDspectra,..This.is true not only for, this simple 
example, -but.for any.X.factored calculation. A 
transition of species A is allowed only if the Mo values 
PieelowOlher species are, constant... Thus .all transitions 
Lot oOsoOr Out of. 2 matrix, SVeGubmererint Clon. members 


of the same subspectrum. 


Because an this matrix,, species A and species X 
have the same number of groups, with the same F values, 
Mersvecies X subspectra,in Fig. 3-3(b) follow.a pattern 
identical. to that of the species A.subspectra just 
described. However, matrices, belonging to the same 
Species A subspectrum never belong. to, the same species X 
subspectrum. This result can be generalized to any number 
Drospecies-. AM y = +1 between matrices, connected by 


species A transitions. Thus these matrices, cannot belong 


to the same subspectrum for any other species. 
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Sige AP Connected Transitions 


The algebraic relationship between the transition 
frequencies of a particular subspectrum, and unknown 
chemical shifts and coupling constants, is useful only 
if’ those transitions can be recognized in the experimental 
Spectrum. Thus subspectral analysis can be applied only 
when the structure of the.observed.spectrum has already 
been determined. Often this requires special experimental 
ks 


techniques, such as the low power double resonance methods 


described below. 


Iso InDOR!° experiment, a reasonably 


In a spin tickling 
isolated transition (referred to below as transition 
iemeterciosen. it as irradiated continuously while a 
second radiofrequency is swept slowly through the remainder 
of the spectrum. The resulting double resonance spectrum 
indicates which other transitions, in the region swept, 


share an energy level with T. All of these transitions 


(rneludins 1) forma set of connected transitions. 


Transitions can be connected in either of the two 
ways illustrated in Fig. 3-4. If the energy level common 
to the two transitions lies between the other two energy 
devels, as in Fig. 3-4(a), the transitions are said to 
be connected progressively; otherwise, as shown in Fig. 
3-4(b), they are connected regressively.!° With careful 


adjustment of experimental conditions, an INDOR spectrum 
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FIG. 3-4. CONNECTED TRANSITION TYPES 
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can show clearly whether a particular transition is connected 
Pro ieesiVely Of Tegressiveryrto- 1. see’ so.. for v-examples. 
Tickling experiments can also give this information in 


Sones Cases’. 


Thus a double resonance spectrum can allow assignment 
of the relative positions of all of the energy levels 
belonging to transitions in the connected set. Further 
double resonance spectra, in which other transitions 
in the set are irradiated can extend the collection of 


energy levels and transitions whose structure is known. 


The two approaches to the structure of a spectrum, 
au terms of subspectra, or in terms of sets of connected 
transitions, are complementary. Subspectral analysis 
leads directly to values for unknown chemical shifts 
and coupling constants, but it requires a knowledge of 
the structure of the spectrum that often can be obtained 
only from double resonance experiments. The pattern 
of connected transitions, in addition to allowing assignment 
to subspectra, often provides direct confirmation of 
the signs of coupling constants in the many cases where 
these signs do not affect the appearance of an ordinary 


Single-resonance spectrum. 


It is seldom possible to work out the entire structure 
of a complicated spectrum from double resonance experiments. 
If other transitions T°, T°~, etc. are coincident in 


frequency with T, every transition connected to any of 
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haben y  ,; ... Will. appear, imatherdouble resonance 
Speettum, to belong to’avsingle connected set. Any 
additional double resonance spectra intended to resolve 
the problem may be equally ambiguous. One can, however, 
use whatever partial structure has already been determined, 
to assign transitions to subspectra. One then calculates 
a trial spectrum from parameters obtained by subspectral 
analysis. If the double resonance results were correctly 
interpreted, the experimental partial structure will fit 
into the calculated complete structure. The calculated 
structure will also suggest useful transitions T for 


additional unambiguous double resonance experiments. 


The relationship between subspectra and sets of 
connected transitions is discussed again in 84.1 in 
connection with the design of a computer program for 
spectral calculations. The analysis of the spectrum 
Ofecissperfiluoro-Z-butene,(ints’5.1;¢illustrates™the use 
of INDOR to determine part of the structure of a complicated 


spectrum. 
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CHAPTER 4. COMPUTER CALCULATION OF SPECTRA 
avtT; OH TeerMve SALor ea Computer Program 


The usual goal in analyzing a high resolution 
NMRe Spectrum” 1s”>to sobtain -all"of-the chéemical@shiftts 
and coupling constants. Computer programs, generally 
based on either LAOCOON!” or NMREN/NMRITL®8 , have become 
an essentaalyaadeto;the.analysis.of complicated. spectra. 
Some of the limitations of existing programs can be avoided 
if a new. program's design is derived from its intended 
uses, which vary according to the type of spectrum being 


analyzed. 


The spectra of many molecules without frame symmetry 
contain several nearly first order multiplets well separated 
from each other. Spacings in these multiplets give 
reasonable approximations to the magnitudes of most or all 
ofthe coupling constants. The-centre’ of gravity of each 
first order multiplet gives an approximate chemical 
Brite. Frequently, particularly in™proton spectra, ~ there 
remain some chemical shifts close enough together that 
only rough estimates of their values can be obtained 
by inspection. To get more accurate values for these 
chemical shifts, and any remaining coupling constants, 

a series of trial spectra are calculated for various 


reasonable values of the parameters. 
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mhe essence Of this type of analysis is obtaining 
a visual match to the observed spectrum. A good plotting 
routine 1S,essential; other characteristics of the program 
become important only if they limit the size of problem 
that can be handled. When there are highly complex over- 
Pappineemultiplets, a plot that simulates experimental 
JanesShapes 1S necessary. Otherwise, a simple stick (bar) 


plot is adequate. 


Visual matching of the observed and calculated 
spectra does not guarantee that the correct parameters 
have been found. A seemingly straightforward multiplet 
structure may be a deceptively simple limiting case of 
a more complicated type of spectrum, ?° Even when the 
Spectitim 15, aS )simple.as.itrappears,ito.be,,.a first, order 
analysis. does not, reveal the,relative. signs, of coupling 
constants. These uncertainties can be resolved only by 
comparing the structure, as well as the appearance, of 
the observed and calculated spectra. This places 
heavier demands on a computer program than does trial- 
and-error visual matching. Thus the program features 
discussed below in connection with the analysis of 


complicated spectra are also useful in the full analysis 


of simple spectra. 


Molecules whose NMR spectral parameters cannot 
be obtained by inspection fall into three classes. The 


first class includes the many common small spin systems, 
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such as (hd and [AX], , for which there are algebraic 


solutions, tables of examples, or both. 1)? 


The, second 
class includes those larger spin systems whose spectra 
are simplified by some combination of symmetry and 

X factoring. Their complicated complete spectra contain 


simpler subspectrat! 


Sitiitar git tOrmato the spectra we 
class 1 molecules; ‘from these subspectra, the chemical 


shifts and coupling constants are calculated. 


The third class, large molecules lacking both 
Syimetive and) X tactoring , is of little concern heresas 
it seldom yields spectra that can be analyzed. The 
positions of the strong, lines do not correspond to .chemical 
efgtts and the Spacings (few ot which are repeated) 
bear no simple relationship to the coupling constants. 
Double resonance experiments (see 83.3) may be of little 
value, as even the sharper lines are very often super- 


positions of several transitions. 


Molecules in the second class, whose spectra can 
be analyzed, but only with difficulty, are the ones which 
most require a carefully designed computer program. 
Usually the most difficult stage in the analysis is the 
assignment of the observed lines to the various subspectra. 
Since the simpler subspectra usually originate in matrices, 


with high statistical weight, some of their lines may 
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be intense enough to stand out clearly. Others will 

remain hidden among the lines belonging to other subspectra. 
The first step in the assignment is often the calculation 
of trial spectra for reasonable starting values of the 
parameters. If a trial spectrum bears any resemblance to 
the real spectrum, it will show which prominent lines in 


the real spectrum belong to useful subspectra. 


To be useful in such subspectral analysis, a program 
must print enough information to allow the convenient 
tracing of calculated subspectra-in two directions. 
Starting from a transition, one needs the quantum numbers 
that define its subspectrum; starting from a‘ set of 
quantum numbers, one needs all of the transitions belonging 
Pomciatesupspectrum. ~The aecont requirement cannot be 
conveniently met by an ordinary transition table arranged 
by frequency, so the program must provide a second table 


arranged by subspectra. 


The weaker transitions in the simple subspectra 
usually cannot be found without using double resonance 
feenodcemassdescribped in@33.3.)-)lf one of the strong lines 
already assigned is sufficiently isolated from other 
lines, it should be chosen as the irradiated transition. 
INDOR can be used to particular advantage, as it reveals 
the frequencies of connected transitions even when they 
are weak components of strong, unresolved multiplets 


teen 85.1). 
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To aid the planning of double resonance experiments, 
and the interpretation of the results, a computer program 
must show the structure of the calculated spectrum in 
Peruseo, Sets of iconnected "transitions. 2 This ‘can. best 
be done by printing lists showing what transitions are 
connected to selected transitions chosen by the user. 

A table giving the transitions connected to every transition 


in the spectrum would be too large for convenient use. 


It has been implicitly assumed above that the 
structure of the calculated spectrum is the same as that 
of the observed one. In a sense this is always true 
if correct parameters have been used in the calculation, 
but unless the program uses all of the factoring appropriate 
Poweach spin system, the structure of the calculated 
Spectrum is not apparent to the user of the program. 
Since a program that ignores valid factoring does not 
have available all of the quantum numbers that define 
the observed subspectra, it provides little aid in a 


detailed analysis of the spectrum. 


A computer program should be designed for the 
convenience of the user, even if this requires more work 
on the part of the programmer. The input data should 
be arranged logically, using a minimum number of different 
data card formats. Errors in the data should produce 
clear messages rather than a cryptic failure of the program. 


There should be a wide variety of output options to allow 
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for the different types of analysis described above. 


Observed and calculated spectra are more readily 
compared if the chemical shifts of each species can be 
given as displacements from the same reference frequency 
used experimentally. Each species' spectrum must then 
be given separately, to avoid confusion between numerically 
Similar frequencies belonging to different species. This 
has the further advantage that each species' calculated 
Soectrum 1a5.readily limited to a particular region: of 


amterest, OT Suppressed entirely. 


Finally, a program should be capable of making 
calculations on large spin systems, and should be 
efficient enough to make such calculations economically 


feasible. 
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At2* The Computer 


The design of any computer program must take into 
account the limitations of generally available computer 
hardware. The most serious limitation is that a computer's 
fast main memory, in which all calculations take place, 
is small compared with the data storage requirements 
of an NMR program for large spin systems. Large calculations 
Caieove Made only if, throughout the program, anything 
Poemrelevant. to. the current.step of the calculation is 
moved out into slower auxilliary storage. Most computers 
allow program instructions as well as data to be moved 


in and out of fast memory as needed. 


An efficient NMR program thus consists of a series 
Boesteps,.one of which is, active in main memory at a 
ime. Each step processes a string of data blocks, typicaldy 
of matrix, SiZ2e in assemo Ly -iinertasitol.o glesnaturds 
Setmor steps, which closely parallels the stages of a 


Pamecatcuiativon.. 15 Siven an (4.35.1. 


There are two main kinds of auxilliary storage. 
Nearly all computers are equipped with sequential auxilliary 
storage, usually magnetic tape. In order to reach a 
particular piece of data, a sequential storage device 
mist eomnMalklwystant, at the beginning and read through 
all of the preceding data. Thus data items can be retrieved, 


conveniently, only in the same order in which they were 
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stored. 


PiLeCt access auxildivary, storage, stich as magnetic 
disk, has become increasingly common. Data items can 
be retrieved in any order, but at the price of increased 
overhead time spent keeping track of where each piece 
of data has been stored. The overhead is reduced if 
a direct access device is divided (logically, not physically) 
into a number of regions. Any region can be chosen at 
random, but within a region only sequential access is 
allowed. This gives the programmer, in effect, dozens 


of fast sequential storage devices. 


An NMR program with full factoring is sufficiently 
large and complicated that it is most often run on a 
large computer in the same general category as the 
IBM 360/67 at the University of Alberta. It is frequently 
more economical, with such a computer, to design a program 
to use a minimum of computing time, even if it then requires 
more main memory space. However, it is also desirable 
for an NMR program to be suitable for use with a small, 
fast computer of the kind now frequently part of an NMR 
spectrometer installation. The program described in 
64.35 Can®be adapted®* to either’kind'of computer, but requires 
Pereeteaccess auxilliary storage if it ds to run efficiently 


on a small computer. 
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Weds Some Aspects of Program Design 
Ago el. introduction 


The NUMAR™ computer program, whose listing is 
given in Appendix B, has been designed to meet most of 
the requirements of 84.1. The present section provides 
a critical discussion of the methods of calculation used 


in NUMAR, and describes some alternative methods. 


The basic structure of NUMAR is determined chiefly 
Byecne use Of all possible factoring, to the matrix, 
level, throughout the program. This makes the program 
more complicated, but saves both main memory space and 
calculation time. The design has also been influenced 


bya decision to use only sequential -auxilliary storage, 


making the program suitable for a wide variety of computers. 


Although they may be combined in various ways, 
essentially the same set of program steps appears in all 
programs that calculate NMR spectra. Arranged in the 
order in which they appear in NUMAR, they are: 


Read and print the input data 

Construct the basis functions 

Generate the Hamiltonian matrix 

Diagonalize the Hamiltonian matrix 

Calculate the transition frequencies and intensities 


Order the transitions by frequency 
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Ono iat the transttron tist 
h Plot the spectrum (bar plot) 


4 Plot the spectrum (full lineshape simulation) 
In the overall flow-chart for NUMAR, Fig. 4-1, 
the blocks are labelled according to this list, as well 
as by the names of the FORTRAN subroutines (see Appendix B). 
The use of full factoring has made most of these steps 
large enough to remain separate, but step c has been 


combined with d, and 4 with g. 


The step whose data storage requirements limit 
puobelenesize is TRANS. It must.deal with a pair of matrices,, 


compared with one matrix, for HAMIL. BASIS deals only 


4 
with vectors of quantum numbers, rather than matrices, 

so its data storage requirements increase much more slowly 
Weeheproblem size. SPOUT also deals: with vectors (of 
PeanmettLon frequencies, etc.) ;“1ts storage requirements 
anerease only linearly with the number of transitions, 

but the number of transitions increases rapidly with 
Baroremesize, For very large, highly factored, spin 
Systems such as the 2Z,5-dimethyl-3S-hexenes (85.2), a 
species' transition list may require much more storage 
than that used for the matrices in TRANS. MHowever this 
need not limit problem size; the transition list can be 


broken up into frequency regions before sorting without 


meauasticyioss of efficiency. 


PUGS 4-1 also shows an additional step, labelled 


ADJUST, for the iterative adjustment of chemical shifts 
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iets NoOcupart Of ‘the currentJNUMAR™ program, but: could 
be added with few changes in the other steps. Possible 


designs for ADJUST are discussed in Appendix E. 
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FIRST CASE 
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SPECTES 
SPECTRUM 


<> NEXT CASE 


FIG, 4-1, NUMAR FLOW CHART, 
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AS. 2 Basis Functions 


The choice of basis for an NMR program determines 
Pie uceree, Of.factoring of the Hamiltonian matrix, and 
consequently the kinds and sizes of spin systems that 
the program can handle. The ‘basis functions' generated 
in the BASIS step of the NUMAR program are sets of quantum 
bulibers that represent product functions of the form 
given by (2.24). In order to provide all three stages 
femcacLOring Outlined in Chapter 3, these basis. functions 
are grouped by symmetry into matrices,, and within each 


matrix, by M and Mg into matrices. 


POeruUuLyYsdetine a -basisstunction of =the vrormepiven 
Pye2. 24), one must have the total spins F of all of 
the magnetically equivalent groups, in addition to the 
A and M values of all invariant groups and symmetrically 
equivalent pairs. For each group or pair each of these 
quantum numbers can take on any of a range of allowed 
values, in steps of 1. Each quantum number's range of 
allowed values is determined by the current value of 
another quantum number. This natural hierarchy of quantum 
numbers leads to a very straightforward procedure for 


generating all of the product functions. 


All possible sets of valid F values for the 
magnetically equivalent groups are generated. For each 


magnetically equivalent group of n nuclei of spin I, 


67 


and matrices, respectively. 


- ‘, 


= = atop dong subarea oo ¥ 


eonigressh mergord. AMM ts WOR 2 


© 


bas -kivton wained Limelt oft To ete 


{3 pineteye wige Fo +ot te bres ehecta: ob 


hetinereg. ‘enobtonn? 2! i oat areca . —_ ts 


mises sbo 2tee vte tint ig SLAB ii anit to: sda o 


ago? sae akesd »acdé .f 1rotgndd 08 bemeteae 


fses aititiw bos ,-2072TIRN OSE YIOSIS 
. ‘. = oat) 7 — 
levity eqest. ,2o3i12 s ,»wooliden egnt gm eee 
> ) 7 
: , a wie ee 
‘ 
ry} Gg 4 5 ‘ E 
:  ) 
‘ + 
CT. oe 
rf ‘as S Lob 
—_ 
Sf i ; zijT tes ft 
‘ 
WED) rie - Oat hi ; fie { ; 
' 
’ ’ : Pit = 
wu ee) yf } ( t > o Sep 


bovolis to sahse sito _ VP Ne ipaédp 
. 7 X ” 


to Q2h61 2*1edmya tmjtabup dost of) ie eapeee 
: F a : = —a 

‘ P 4 e yo. 

to Sofsv Jnetye> sn? ¥r- leah ae ee as hE. ; 

4 sees - 

‘ vada x Wp 
. i . a a ih 

T6462 oeuboto1g brewro? ai 8 EG ied WE OM f 
— 


—ait yo} aotfsy Tohilev 


foes 07 bo lax og ors =equatg 


el digs te fefoum a Io 


the range of valid F values is given by 

od ee ee OTs LE = ni. (451) 
For each of the sets of F values all possible combinations 
of A values for the invariant groups and symmetrically 
equivalent pairs are generated, subject to the constraints 


Bee geese Est | (4.2a) 


fOr cach pair, and 


> 
Ht 
Ls | 


(4.2b) 


Doaweech Invariant group. “For each- of. the’ sets’ ofA 
values, all possible combinations of M values are generated, 


Subject to the constraints 


A2M2 -A (4.3) 


for each pair or invariant group. 


It is convenient to generate the basis functions 
for one matrix, at a timess All on ‘thesbasis tunct1ons 
penerated for a particular-set of F values belong to 


one magnetic equivalence submatrix. If the molecule 


has no frame symmetry, or if F # F° for any pair, they 
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are the basis functions for one matrix,. 


When the molecule has twofold frame symmetry and 
F = F~ for all pairs, the magnetic equivalence submatrix 
is factored into symmetric and antisymmetric matrices,. 
The basis functions generated according to (4.3) for a 
particular set of A values have a symmetry given by the 
paoduet ‘overvalivpaits 


AE Cr aA (LF (4.4) 


Cs) 

fg 
When basis functions in a Symmetric matrix, are being 
generated, the antisymmetric sets of A values generated 
according to (4.2) are skipped; when the basis functions 


mrethe corresponding antisymmetric matrix. are being 


Z 
MPenerated, the sets of A. values just used for the symmetric 


matrix, Poees kK LOped. 


Successive basis functions produced by this generating 
procedure seldom belong to the same matrix, or even the 
same matrix,. They must therefore be sorted by M into 
matrices , and awd thin, each matrix, bY <a11L. of the Mo 
values into matrices,. To, ald this «sorting procedure, 
M and Mg values as well as A and M values are included 


as part of each basis function. 


In the program, each of these spin quantum numbers 
is represented by an integer equal to twice its true 


value. Memory space is saved by 'packing' the integer 
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Fepresentations of. ‘the Mg values of each basis function 
into a single integer, and those of the A values into 
another (see Appendix C). Because each possible set 
of My values within a matrix, has a unique packed represen- 
tation, one need only check whether the packed Mg values 
are equal in order to determine whether two basis functions 
belong to the same matYix,. 
This allows grouping of the basis functions by 
matrices. and matrices, using a simple single-stage sorting 
routine. The basis functions are ordered by decreasing 
M value or, when M values are equal, by decreasing packed 


Mo value. 


The sorted basis functions for the matrix, are 
Pren placed ™@in auxilTtary "storage =" For “each matrix,, 
the F values and the statistical weight are stored; for 
each matrix., the “M-value;> for “each matrix, , the packed 
Mo value. The quantum numbers stored for each basis 
function within a matrix, are the> packed A.value, and 
the M values of all groups and pairs. These quantum 


numbers, forall matrices,, ave ssUtL IT CLen te rosdaet ime 


all following steps in the calculation. 


The procedure used in the ean: NUMAR program 
differs in one major respect from the one just described. 
To simplify the calculation of some Hamiltonian matrix 
elements in those matrices, where Ee =P for allepairs 


(see HAMIL, §4.3.3), NUMAR reverts to a simple magnetic 
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equivalence basis in the other matrices. where frame 


symmetry provides no additional factoring. 


Heit IN put tO NUMAR each pair is tredted as 
a single entity. This representation is then expanded 
internally (see Fig. 4-2) to give two magnetically equivalent 
groups in place of each pair. When F # F~ for any pair, 
these two groups are treated in the same way as any other 
Maonetically equivalent groups. ,-When F = Fe for all 
Paats, the first yf the/two groups (for each pair) is 
given the characteristics of the pair, and the second 
@segiven spin /F =/0 in order to effectively remove it 


from the calculation. 


The dual numbering system makes NUMAR considerably 
more difficult to understand and creates complications 
in HAMIL nearly as great as those it eliminates. The 
simpler approach described first, which always treats 
a pair as a single entity, is recommended for future 


programs. 
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eh. O's The Hamiltonian Matrix 

The overall design of this step in the calculation 
1s straightforward. The basis functions for one matrix, 
are brought into main memory, the matrix elements are 
calculated, the matrix is diagonalized, and the eigenvalues 
and eigenvectors are placed in auxilliary storage. When 
there is frame symmetry, some details of the calculation 
vary according to the procedure chosen for handling those 
matrices, that cannot be factored by frame symmetry. 
The recommended procedure is described first, followed 


by a brief summary of the ways in which the NUMAR procedure 


Ppreters  fromsit. 


All Hamiltonian matrix elements are calculated 
using the Hamiltonian (2.23), which requires that the 
symmetrized combinations of coupling constants given by 
(2.22) be used in place of the actual coupling constants 


connecting symmetrically equivalent pairs. It is assumed 


RS: 


that a coupling constant matrix containing these symmetrized 


coupling.constants .has.been;createdpin;ythe pinputystep, 

In this-matrix, the. only nonzero elements, below the main 
diagonal are the difference couplings Ine between pairs, 
and the only nonzero elements on the main diagonal are 

the couplings Jop- between the two magnetically equivalent 
groups in a pair. Each symmetrically equivalent pair 


isetreated, at all times, as.a single entity; the size 
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of the coupling constant matrix is given by the number 


of pairs plus the number of invariant groups. 


Before any matrix elements are calculated in each 
matrix,, its F values are brought into main memory. Two 
vectors are needed, the second containing only the F~ 
values for the symmetrically equivalent pairs. Most 
matrix elements are calculated using only the A and M 
values of the individual basis functions (A = F only 
for invariant groups), but F and F* are used directly 
in calculating the matrix elements of ReR- and B 


for the symmetrically equivalent pairs (see 82.3.2). 


Rather than waste time evaluating FR, and the 
Seoressions for 4(F,F ,A) and g(F,F°,A) (see Eq. 2-16) 
every time they pappear Ain Jartmatrix ielement, fone icantcaliculate 
them, for each pair, ati the beginning of each matrix, 
and store the values for use when needed. The values 
On 40h 2A) pand tq (F pF7,A) afexstoted incmatrices.; 
indexed by pair and by A value; the values of F+R? are 
stored in a vector indexed by pair. The rows, in the 
4 and g matrices, corresponding to invariant groups can 


beset to zero, since they,will never be used. 


basis ptunctions for éa matrix, are brought into 
main memory; each basis function consists of integers 
representing the M values of all pairs and invariant 
groups, and a single integer containing the corresponding 


A values in packed form. The Hamiltonian matrix elements 
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Hae are then calculated, one row at a time. 


The first step in calculating the matrix elements 
of row i is the unpacking of the A values of basis function 
i. The diagonal matrix element Hes is then calculated; 
each diagonal element is a sum of terms contributed by 


meeieparts 01. (2.23) . 


since the Hamiltonian matrix is real and hermitian 
(4.€. Symmetric), one need calculate only those off-diagonal 
matrix elements He that lie above the main diagonal. 
tiie tne last two sums in (2.23) can provide off-diagonal 
matrix elements. If the packed A values for basis functions 
i and j are identical, any nonzero off-diagonal element 


1s a member of the A sum; otherwise it is from the B sum, 


When the packed A; and ae Values taresnot tdentical., 
the A values for basis function j must be unpacked. His 
Pewee ro unless the conditions: (2.14) ‘are met for two-pairs, 
for example, pairs p and q. Combining these conditions 


with others inherent in the factoring into matrices, 


and matrices,, one requires for a nonzero matrix element 
erin -B : 
qp2ip)*B(a) 
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and j.) If basis functions i and j meet these conditions, 
the matrix element Le 1s calculated using (2.15) and 


the stored values of he seep ced pe 


When the packed A values are equal, there is no 
need to unpack the A values of basis function iach vel an 
conditions on AM similar to those above are satisfied, 

a matrix element AG = J yqh lp) Ata) is calculated using 
(2.13). Note that all nonzero off-diagonal matrix elements 
involving invariant groups (including those where p is 
@eiivartant group, dand’'q'* “is ‘a™pair) belong to the A 


sum, not to the B SUEY jee te earn 5 ee 


When the matrix elements in all rows of a matrix, 
have been calculated, the matrix is diagonalized. The 
diagonalization routine used in NUMAR, based on a new 
variant of Householder's merhedte.t is particularly efficient 
in its use of memory space, very fast, and very accurate. 

It has already been pointed out that at least two matrices, 
must be present in main memory during the transition 

Par culLlation. lhus the use ‘ot ‘double precision arithmetic 

in the diagonalization does not increase the total amount 

of data storage used by the program. It has the advantage 
that the single precision eigenvalues and eigenvectors 


placed in auxilliary storage contain no rounding errors 


from the diagonalization. 


The procedure used in the HAMIL subroutine of 


NUMAR for generating Hamiltonian matrix elements is generally 
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PueravecOscthe one just described. | Akl of the differences 
are the direct result of a decision to use the frame symmetry 
basis (2.24) only in those matrices, factored by frame 


symmetry. 


In the symmetry factored matrices telaAteo Tx el 


Wet 
Sylmecrized coupling constants, Similar tO the one described 


above, is needed. However, in the other matrices the 


9? 
coupling constants must be restored to a non-symmetrized 
form. In HAMIL, before any matrix elements of a matrix, 
are calculated, its«coupling constant matrix is generated 
from a simple matrix of coupling constants created by 
miemriput routine. The generating Iprocess is complicated 
by NUMAR's expanded internal representation of a pair 

as two magnetically equivalent groups, which increases 
the size of the coupling constant matrix to the number 
Oreinvariant groups plus*twice=the-number of pairs. The 


relationships among the three forms of the coupling constant 


Poacyixesare, 1 llustrated jin Fig. 4-3. 


In NUMAR, there is no need to generate values of 
g(F,F~,A), since they would all be zero in those matrices. 
for which all F = F~, and symmetry factoring is used. 
Tteigualso unnecessary.ito génerate {(F,F° ,A) for each 
pair at the beginning of those matrices, in which it 
is used. A single matrix of values of 4(F,F~,A), indexed 
by F (=F*) and A, is used for all groups in those matrices, 


where it is needed. The values are calculated at the 
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Depinning Of the program using (2.17). 


Because g(F,F°,A) is zero, no provision is made 
in NUMAR for terms from the B SUN -OL@.( Anco) 21) ae Olas 
matrix elements He. OtHerwrse.. the acttiad process. UL 
generating the Hamiltonian matrix elements is as described 
Poove. In the matrices, that are symmetry factored, the 
second position assigned to each symmetrically equivalent 
pair must be ignored. Because each such position has 
Peeneset to M= A =F = 0, this requires no special action. 
Mion NoOespoecial (procedure Wseneededi for matrices, that 
are not symmetry factored. Since each group always has 
Peer, all the basis. functions an each matrix, have! identical 
packed A values. Thus the procedures’ for ‘generating 
matrix elements from the B sum of (2.23) are never used, 
and the calculation reduces to the usual magnetic equivalence 
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Baus, As. Transitions 


Transition calculations in NUMAR include essentially 
the same steps used in other programs, 1/»18 However, 
X factoring creates some complications in the use of 
auxilliary storage devices. Three approaches are discussed 
below: one requiring direct-access storage, one for 
small computers without direct-access storage, and the 


one used in NUMAR. 


Mathematically, the transition calculation is 
straightforward. The matrix elements of the transition 
PpetetoryxX of Eq... (2.21)s,arexacalculated in thesbasis 
representation, using the selection rules discussed in 
coer elhesXematrix is thenetransfonmedttoathe represen- 
tation in which the Hamiltonian matrix is diagonal, using 
the eigenvector matrix U.obtained previously in the 
meaconajization step. Since U is-real and unitary, its 
transpose eas Squall ecoOL1tSe inverse ie The unweighted 
intensity Ts; OLethestransi tion connecting senereyeleveis 
i and j is given by the square of the transformed matrix 


element: 


_ 2 


This intensity must be multiplied by the statistical 


weight W.. appropriate to the matrix, containing energy 
1j Z 


80 


v2 . : 
¢ 
© P an -* 7 
: a 
: fine 
(iieitesaes obulont AAMDY om: oo 1 EAFOER ee colt bas 
. “; : a Rac 


4 


T sy oor a smergotq talto nb 
Yo. cen ofly. wi. anorteo ii qmes ric 
meyer th oxs aoiisso ng n'y it 4 ) ey" ee a9 ore et x 

bid wee snisiwer aap 
sodttot2 gegote-t2ebEe swale tw wey #09 


. _ 
rd ~ ra 
; 
* = 
* 
7 ‘7 F ? y) Ff rit f ' et Linabsa 
- i 4 a? ‘ ) 
* ge 
_— m | m +, a ¢ : toy 6 af 't 
(moe Tih |e t ets Yo ettimolse xinowey eee 
= 
} vii oko $ 
i bs i ‘ ri Bae \ 
aT ed : (?- ca 


sofeu , fafibpeit ei; ceri t ews fags Dime ody, 


* “4 oe te a 
tid A rr toe Pe 
“ = sa ~~ 


.. 
Po 


2+] (thtinu bite fsa, ebyy, pues ; 
se 2 Cae | ? q 4 “ie = 
bor phewsu ent §.°  “e21 255° ot. pane 


") eee 
elayvet Yorens anidoannme5 nodecanead, ont , va 


cintam berrrotensirt oF to orsupe, ails véon ee ae 


ee nae | oe 
‘ | septs ioees silt 4d botkgte tiny 


| atone ‘yitatetno> pee | 
| ie et Alea 
a a koe 


i 


81 


levels i and j. The corresponding frequency v.45 is simply 


the algebraic difference in energy of the two eigenvalues: 
Vieg=eeE. ©) E« (4.6) 


The NUMAR program follows the usual convention? of calculating 


the emission (AM = -1) spectrum. 


The nonzero elements of the X matrix are in biocks 
that do not fall along the main diagonal, but rather 
connect basis functions in matrices , that obey the selection 
Pulesringss.Z.ceBecausestheyX matrixeisrsymmetric,! only 
the blocks above the main diagonal need be considered. 
The transformation matrix U is block-diagonal, with each 
block corresponding to a Hamiltonian matrix,. Thus the 
transformation of a rectangular block Xnq ofexXsiconnecting 
matrices, Pend geocan.berxwritten: 
vp “pq “a ie 
Note that although Xoq is normally rectangular, as much 
memory space must be allotted to it as to the U matrices, 


because either matrix, p or matrix, q. may be the larger. 


This series of matrix multiplications is the critical 
step in limiting the size of problem that can be handled 
by an NMR program. Part of the available data storage 


space is used to hold a variety of small items such as 
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the eigenvalues of the two matrices, Dp taitd, qie Even” wsaang 

a minimum-storage method of matrix multiplication, the 
transformation of Xd requires space for one other matrix 
of equal size. Thus the program can handle only those 

spin systems whose largest matrix, accupi es less than 

half of the data storage space available for the transition 


Galculation. 


A transition calculation without X factoring is 
straightforward. For each matrix, , the eigenvalues and 
ik vector a i i = = . 
eigenvectors of the matrices. with M Mn and M Me ots 1 


are brought into main memory, and the transitions from 


M to M 
m 


fase aa are calculated, as described above. The 


a 
eigenvalues and eigenvectors for M = 1 ae are no longer 
needed, so their region of main memory is reused for the 
eigenvalues and eigenvectors of the matrix, with M = Meax 2° 


The transitions from M -1 to M —Pitave Utheneca ucilated, 
max max 


a 
The procedure is repeated, with eigenvalues and eigenvectors 
of each matrix, Ssernving stirsteas the ‘lower then the 
Dipper‘ set. Each set.is brought in from auxilliary 
storage only once, in the same order of decreasing M 


in which the matrices, were created and diagonalized in 


the previous step. 


Referring to Fig. 3-3(a), one can see that in 
eaexe factored transition calculation, a,procedure similar 
to the one just described can be used to calculate each 


complete species A subspectrum in turn. As each matrix, 
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belongs to only one species A subspectrum, the only compli- 
cation introduced is the need to bring the eigenvalues 
and eigenvectors into main storage in order of decreasing 


Ma» for each subspectrum in turn. 


If only species A transitions are of interest, 
one can generate and diagonalize the Hamiltonian matrices, 
in the order needed for the transition calculation, and 
retrieve them from sequential auxilliary storage without 
difficulty. However, a completely different order of 
matrices, isqrequired for calculations on éach additional 


species (for example, see Fig. 3-3(b) ). Three solutions 


to this problem are given below. 


When direct access auxilliary storage is available, 
one can generate and diagonalize the Hamiltonian matrices, 
in any convenient order, then retrieve them in the order 
required for transition calculations on each species. 
Tables are kept showing the locations in auxilliary storage 
of the basis functions, eigenvalues, and eigenvectors 
that correspond to each set of Mo values. Ties east avics 
can themselves be placed in auxilliary storage, so that 


Main memory space need be reserved only for the entries 


for one matYrix,. 


The majority of the data area in main memory is 
used for two regions each the size of the largest matrix, 
that the program can handle; they are referred to below 


as regions one and two. 
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Basis functions for the current pair of matrices,, 
p and q, are brought into region one, and used to construct 
the matrix Xd (see Eq. (4.7) ) in region two. The 
eigenvalues of the two matrices, are also brought into 
main memory. The eigenvector matrix U. is then brought 

E ~ 
into region one, and the matrix product P = Unxng is 
formed as follows. The first column of x is used to 
calculate the first column of P, which is stored as a 
vector. This P vector is next copied into the first 
column of Xnq? and then is treated as the second column 
eeeeeeac Ciuc end ofthe matrix multiplication, Xoq has 
been completely replaced by P. The second eigenvector 


matrix Ug is then brought into auxilliary storage, and 


used to calculate Pua 


When each element of Ae has™been™calcutateds “1t 
is squared and multiplied by the statistical weight W 
Poeeive 7a transition intensity. “~“The* frequency 1s calculated 
from the difference between the corresponding eigenvalues, 
and placed in auxilliary storage along with the intensity, 
subspectral quantum numbers, and an origin number to 
indicate which energy levels the transition connects. 
mmiemconvenient to edit.the list _of -calculated transitions, 
placing in auxilliary storage only those whose intensity 
exceeds a minimum value specified as an input parameter. 


Frequency limits may also be imposed. 
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Aeprocedure similar to the one just described 
can be used to calculate any one species' transitions 
when only sequential auxilliary storage is available. 
The basis functions, eigenvalues, and eigenvectors are 
placed in auxilliary storage in the same order in which 
Peeyearcenceded in the transition calculation. {If 4 
matrix, Participates as both ‘lower’ and ‘upper® in the 
transition calculation, its eigenvectors are placed in 
auxilliary storage twice in succession.) When all 
PransitaOns Of one species Have been calculated, the 


entire program is repeated for each additional species. 


Poemthirndsmethod of calculation, used in NUMAR, 
is the fastest of the three, as it calculates the eigen- 
values and eigenvectors of each matrix, only once, and 
brings them back into main memory only once. It does 
Porcerecuire direct access. auxillidry storage, but it 
does require far more main memory space than either of 


the first two methods. 


The basis functions, eigenvalues and eigenvectors 
Gieed Mito f tthe matrices, inftdn-adj dcentipaur of matrices, 
are brought into main memory. A pair of matrices,, p and q, 
Connected by transitions of the first species, is selected, 
and the matrix Xoq is sconstructed,  .Une row of thewprocuce 
matrix P = UXnq is calculated, and Storedfas a vector. 
This row of P is used immediately to calculate, one at 


a time, the elements of the corresponding row of the 
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product ek These are used, as above, to calculate 
transition intensities, and the transitions of significant 
intensity are placed in auxilliary storage. The next 
row of P is then calculated (and stored in the same vector 


Meed. nor the first row), and so on. 


Pinet rds tliOns ioLeall=speciessare calculated tor 
Pris pair OL matrices. before moving on to the text paix. 
Because the eigenvectors etc. of each matrix, are used 
again for each additional species, their storage space 
cannot be reused for another purpose. .Thus space must 
be available for all.of the eigenvectors, eigenvalues, 
and basis functions of two matrices,, and for the matrix 
Xoq" This is far more storage than is needed for the 
first two methods, but less than would be needed without 
feraccoring. Taking the matrices, with M = 1 and M = 0 
Seerie. 5-1 as an example, one would need 72 storage 
locations for the two 6x6 eigenvector matrices U, and 
Ua if there were no X factoring. The eigenvectors for 
the six matrices, in these matrices, require only 32 


Storage locations. In addition, the matrix xq is the 


Size ofa matrix, rather than 2 matrix,. 


In the first two methods, all transitions of one 
species are calculated before any transitions of the 
next species. Within a species, all transitions of one 
subspectrum are calculated before any transitions of the 


next subspectrum. Thus one can readily print the transitions 
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by species and by subspectrum in the output step. In 
the third method, the transitions must be sorted in a 
later step for this purpose. However, the third method 
is better suited to grouping together transitions of 


ee speciesmconnectedyto«a commonftransition; 


Uretnie "three *methods of *transttron calculation 
discussed, the first clearly offers the best compromise 
between speed and efficient use of main memory, and should 
be used for those computers that have the necessary direct 
access storage. The second method should be used only 
when limitations of computer hardware make both of the 


other two impossible. 
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avo. 5S! Input and Output 


The input data cards! 


for NUMAR are designed to 

place logically related data items on the same card, 

and to minimize the number of different kinds of data 

card required. Where possible, the most important items 
erepulacedsat the beginning iof, a card, sand reasonable 
default values are provided if the remainder of the card 

ts lett blank. If an error is+found in -the.data,.the 
program prints a message and skips to the data for the next 


melecule., Chiefly because of the .error,checking, the 


input step is large enough to be kept separate. 


The NUMAR output, much of which is optional, is 
patmtenwethroughout "the *program]* The INPUTSsteép “(Fig 4-2) 
prints chemical shifts, coupling constants, etc. on which 
the calculation is based. BASIS prints the total transition 
intensity, and the maximum possible number of transitions, 
Eorpeaci=species. “HAMIL prints the energy levels, including 
in the table all of the quantum numbers necessary to 
define each matrix, , matrix,, and matrix,. TRANS prints 
a Key, to show which energy levels are connected by each 
mene ition..... SPOUT. prints, .a transition, list for each 


species, including all of the subspectral quantum numbers 


: The format of the data cards is fully described in 


Appendix D. 
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for each transition. Transitions are ordered by transition 
number in the Key and by frequency in the SPOUT list; 


both number and frequency appear in both lists to simplify 


cross reference. 


This output does not fully meet the requirements 
ee teiyedsethnere 1S nO Provision 1oOr printing transPtLons 
by subspectra, or by connected sets. However, it is quite 
Gasy tO tind all of the transitions in a subspectrum. 
One first finds in the Energy Level table the matrices, 
With the correct set of quantum numbers, then finds in 
Peeeseye the transitions connecting those matrices,. All 
Pransitions Connecting a partictlar-pair of matrices, 
appear in sequence. In contrast, transitions connected 
to a particular transition T are not in sequence; one 
must search a region of the Key, in both directions from 
T, for transitions of any species whose upper or lower 
energy level is the same as either of the energy levels 


connected by 7. 


In NUMAR, Lie 1Tens Lions sOl eat lespecl essa lonp laced 
in a single sequential auxilliary storage unit as they 
are generated. In SPOUT, this entire list must be searched 
once for each species whose transitions are to be printed. 
If enough auxilliary storage units are avalLlabiegpat 
is more efficient to’use a different storage unit for 


each species' transitions. 
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The number of transitions that can be sorted in 
SPOUT for each species is limited by the amount of main 
memory space available. Some large molecules (such as 
the 2,5-dimethyl-3-hexenes, §5.2), have many thousands 
of transitions of significant intensity within a single 
species a;yOne Ganrvedit.the.transition, list in PRANS, 
keeping only the transitions within a limited frequency 
range. This is an excellent solution when only a limited 
frequency range is of interest; it is inefficient when 
it requires repeat calculations to obtain the remainder 
Oietneaspectrum...ilt;would;be; better to calculate.ald 
the transitions of interest, separating them into frequency 


ranges in TRANS if necessary. 


Rebar pilot, produced onthe computer's output 
printer, may be more useful in preliminary calculations 
frzeeither a transition list or’a full lineshape plot, 
because it provides both frequencies and intensities 
for all the well-resolved lines in a convenient form. 

A useful feature of the NUMAR routine, BPLOT, is the 
option of plotting the entire spectrum of a species 
with any large regions of blank baseline automatically 


omitted. 


A plot that simulates experimental lineshapes is 
better when one is attempting to fit a region of overlapping 
multiplets. The CPLOT routine in NUMAR provides a variety 
of lineshapes and scaling options to allow simulation 


of a variety of experimental conditions. 
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Gra lok 5% NMR SPECTRA OF SOME SYMMETRIC MOLECULES 


This Chapter presents analyses of the NMR spectra 

of a number of molecules containing symmetrically equivalent 
pairs of magnetically equivalent groups of nuclei. The 

full analysis of the spectra of the perfluoro-2-butenes 
(85.1) illustrates most of the techniques of spectral 
Boeaiysis described in §3.3 and §4.1.- The 2,5-dimethyl- 
3-hexenes (85.2) and the trifluoromethyl phosphorus compounds 
(85.3) are typical molecules of current interest whose 
spectra are not readily calculated except by using the 


NUMAR computer program (Appendix B). 
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Sed The Perfluoro-2-butenes 
vd ee Introduction 


Fluorine coupling constants have remained a subject 
of interest for some years. A variety of evidence has 
been found? for a 'through-space' coupling mechanism, 
said to be responsible for large coupling constants between 
*°F nuclei in close proximity spatially, but separated 
by many bonds. More recently, liquid crystal work?" has 
shown that some '?F-'°F indirect coupling constants are 


highly anisotropic. 


Cis-perfluoro-2-butene was chosen as a small molecule 
of high symmetry, one of whose coupling constants might 
have a through-space contribution. Its spectrum (in 
an isotropic medium) was fully analyzed with the aid 
of the INDOR double-resonance technique. The spectrum 
of trans -perfluoro-2-butene, for which no through-space 


effects were anticipated, was also analyzed. 
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Be ie Experimental 


A mixture of cis-and trans-perfluoro-2-butenes 
(Prerce Chemical Co.,«Rockford, Ill.) was introduced into 
a vacuum system and degassed by the freeze-pump-thaw 
method until the residual pressure at liquid nitrogen 


temperature was less than two microns. 


Several attempts to separate the mixture by gas 
chromatography failed. Each column listed in Table 5-1 
was used, at various temperatures from 25° C. to the 
maximum listed in the table, with helium carrier gas at 
flow rates from 30 to 300 mil./min. A Gow-Mac TR2-B thermal 
conductivity detector was used throughout. Under some 
conditions, small peaks indicative of trace impurities 
could be seen in the chromatogram, but the major components 


always gave a single, large peak. 


; ‘ 2 
Separation of these isomers has been reported. : 
However, column 3, used under conditions essentially the 


we gave no better 


same as those of Bright and Matula, 
results than had been obtained with columns 1 and 2. 

Bright and Matula did, however, condition their column 

with ‘three 250 torr samples of COF,' beforesusing site 

No attempt was made to learn more details of this procedure, 


in order to repeat it, since the NMR spectra of both 


isomers could be obtained from the mixture. 


All NMR spectra were obtained from a sample consisting 
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of 0.30 ml. of the perfluoro-2-butene mixture and tics Eas wie 
of CFC1., both previously degassed, sealed into a 5 mm. 
O.D. medium wall NMR tube. Liquid volumes were measured 
at room temperature, in a calibrated tube on the vacuum 


rack, before distillation into the NMR tube. 


‘°F spectra used in the analysis were obtained 
at 14.092 kgauss (56.4 MHz.) on a Varian HA60I spectrometer, 
at 21.138 kgauss (84.7 MHz.) on a Bruker HFX-90 spectrometer, 
anuvat 25.487 Kgauss (94.1: MHz.) on a Varian. HA-100-T12 
spectrometer. +°F INDOR spectra were obtained on the 
HA-60I and HFX-90 spectrometers. Modifications to the 


Varian spectrometers are described in Appendix E. 


All spectra were obtained at ambient probe temper- 
@eures) (300°-350° K) at a scale of 1.0 Hz/cm,7 and a sweep 
rate of 0.02 Hz./sec. or slower. The spectra were calibrated 


as described in Appendix F. 


Each observed multiplet is symmetrical about its 
chemical shift. Thus the mean of the frequencies of 
each pair of equivalent peaks on opposite sides of a 
multiplet is the centre frequency of the Wibod OLens 
(Because of the finite sweep rate, this value 15 -slightiy 
different from the true chemical shift.) For each multiplet, 
the RMS deviation of these values (from their mean) was 
less than 0.04 Hz., the larger deviations coming from 
the weaker peaks. This value can be taken as a reasonable 


estimate of the error in the observed frequencies. 
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nee 5% Assignment and Analysis of Spectra 


A preliminary spectrum of the mixture sample revealed 
four multiplets centred at “100c75  ~145.00 5 "69. o, rand 
67.4 ppm from Scien with relative intensities approximately 
53:2:9:6. These intensities immediately suggested the 
assignment =C-F (isomer A), =C-F(B), =C-CF,(A), =C-CF,(B), 
respectively, for the four multiplets. Comparison of 
the chemical shifts with those of the 2-chloroheptafluoro- 
2-butenes 24 Supported this assignment, and allowed tentative 
identification of isomer A as trans- and B as cis-pert luoro- 
Z2-butene. The coupling constants subsequently obtained 
from a full analysis of the spectra are consistent with 


‘this assignment of the multiplets. 


No strong line in any multiplet is more than 40 Hz. 
from the chemical shifte- Atyl4.dikeaussi othe. two most 
elosely,spaced multiplets (those at.-69.8 and. -67.4.ppm.) 
are more than.135_Hz. .apart;,at.the higher,field strengths, 
the spacing is correspondingly greater. Thus the use | 
of a sample containing both isomers does not interfere 


With the analysis. 
i A positive chemical shift indicates a resonance frequency 
higher than that of the reference, at fixed magnetic 

field strength. A negative chemical shift indicates 

a resonance field strength higher than that of the reference, 
at fixed frequency. Spectra are presented with chemical 
shift decreasing algebraically from left to right, unless 


otherwise noted. 
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A few small impurity peaks were noted in a preliminary 
Survey spectrum. One doublet of separation 3.95 Hz. 
1s only 0.195 ppm. to high field of the cis-CF, chemical 
Siitt. Lt 1S weak. (about 1% of the intensity of the 
cis-isomer), and does not interfere with the analysis 
because the lines it obscures at 14.1 kgauss are visible 
at 21.1 kgauss, and vice versa. The other impurity peaks, 


which are equally small, are far from any of the four 


Mmaitiplets of interest. 


Bomecachumu kia plet:, sspectrea obtained eat, L415 21.1, 
and 23.5 kgauss can be superimposed.. This lack of effects 
dependent on relative chemical shift indicates that 


the =C-F and =C-CF, resonances are far enough apart, even 


3 
moe Keauss, to be treated as. separate species. 

X factoring has therefore been used throughout 
the analysis below, with =C-F labelled species 1 and 
=C-CF. labelled species 2. A common criterion for the 
Hoewo. XX. factoring in a.calculation, is that the absolute 
value of the ratio J/A (see 83.1) between species be less 


u The largest |J/A| ratio in’these molecules, in 


than 0.1. 
Premworst case, (at 14.1 *kgauss} Zs ‘0°002-for*the cis- 


isomer, and 0.004 for the trans-isomer. 


The structural formulae of the two isomers suggested 
that each should have twofold frame symmetry. Since 
the three fluorines of a CF, group are normally. magnetically 


equivalent, both isomers were assumed to be [AX], spin 
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systems. Because each observed multiplet is symmetrical 
Bpoue its chemical shift, each’ chemical shift was chosen 
faethe zero of frequency for its species. This left 

only the four possible coupling constants to be determined 


in the analysis of each spectrum. 


These four coupling constants are labelled 
rake manner consistent with the input to the NUMAR computer 
program.(see Appendix D). For each isomer, the coupling 
constant between the two symmetrically equivalent olefinic 
fluorines is labelled Jay: and that between the two 
CF. groups Jo: The three bond coupling constant F-C(=)-CF, 


is Jio> and the four bond coupling constant F-C=C-CF, 


is Joy: It is convenient to use the tee eis 
SS Oey oe ney el 
Cat} 
Se Gaede N = Jyot494 


weight 4, with all composite particle spins Fi= 4%. Its 
subspectra are identical to those of the [AX] , spin 
system, for which an algebraic solution is available 
(pst, §6-7). The magnetic equivalence submatrix is 


factored by frame symmetry into symmetric and antisymmetric 


matrices,, of dimensions 10 and 6 respectively. The 
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symmetric matrix, has been discussed in §3.3; the antisymmetric 


matrix, has three matrices,, aia 2 © Withee eele Uses 
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Mute Aotactcoring, the matrices x with M = +1 are further 
factored to 1x1 matrices), and the matrix, with M = 0 
Sentains a ‘single 2x2 matrix). 
The subspectra belonging to these two matrices, 

are identical in the species A spectrum and the species 

X spectrum. The most prominent feature is a doublet 

of separation N which originates in the symmetric matrix,. 
in the species A spectrum. this doublet may not stand 
Geceamong the. lines from other subspectra. In the species 
wespectrum it 1s part of a doublet’ of separation N that 
contains half of the total species X intensity. This 
doublet, which includes all the lines from subspectra 

with My emer isecharacteristic Of. the ispecies + spect rcum 


oad 1 [AX], spin systems.°> 


No coupling constant combinations other than N 
are available by inspection. However, analysis of the 
spectrum of the trans-isomer was made simple by the 
availability of literature values for coupling constants 
in a wide variety of related spin systems. A least-squares 


Borrelation of substituent effects on Je_c=c-F gave an 


estimate-° Ping = 139-4 Hz. From the spectrum of trans- 
Pechiioroheptafluoro-2-butene-- lJ 541 = lo Hee and 
lJ55| = 5.45 Hz. From various phenyl substituted perfluoro- 


alkenes’ |J,4| etd , SuHERS lJ52! = ¥.toalgnuat, tand 


: 2 
Ww = "93 to 28 Hz. Probable absolute signs are 5 


21 


Pig ws a2zaeite2i7 


ui {) 


‘ 
=> 


(ebeont ee 
vw 3e 


bitnuin owt s2edd of Sn ainbod avtoey 


¢ a ty, Bo ath me belo Ls »£3OHs ody ip 

re | Ah ox 
sailed. nhs aot ues) $righe mh tem pili aw: 
BAL “f Seay SA4 nei eet eart a | oO: bathe all Bs 
. » an 
ut JeTaneh a; mers 29g A ih 
. yy 
ee wp wat 

a4 chet dine , 4S ape weskt of f° 


= 7 "ies 
LHmOs a6? Bhat 
) qd 1 


i (ie zebutoml ‘Wok 


- 
si 


- pie ie =. we $ 7 : | tats y3 5B TEND ak : 


- ts ; ‘tL rtinga eel +2016 2 gal gan 
io 2iavinos. .tevewee) eee ies 


i) 


200 23 ' TI ig r hod cae t3U75¥ scwseenst § 


qosaine-sessl Aare ana. tbage boron ve ‘tend ae 
é i ee 
A5-SV8S. wiacale Pen 2590148 saapn heehee 9 
wae ry iro 
7 ; Pa ee : 7 fog 
“@5bT2 10 Muy Oege eas owt “EH re cs a =\g ft 


+ a's wt 


>a 


Trial spectra were calculated using values from 
the above ranges, consistent with the experimental coupling 
Boustant sum |.N| = 12.80 Hz. A few trial calculations 
with the NUMAR computer program (Appendix B) yielded a 
good visual match to the observed Spectra with parameters 
Jay Sela 7. , J55 eb) Leer Ji mo. 2a? fae 
Joy = +21.9Hz. No good match to the observed spectrum 
could be obtained if Jay and J55 were given the same 


sign. 


These values were used as initial parameters for 
moeiterative calculation with aie LAME! computer program. 
The final parameters in Table 5-2 give an RMS deviation 
of less than 0.05 Hz. between observed and calculated 
frequencies. Spectra calculated from these parameters, 
and the corresponding observed spectra, are shown in 
Preeori. The ottter®lines of the CF region are omitted, 


as no satisfactory experimental spectrum is available. 


An INDOR experiment which confirmed the relative 
Signs of the coupling constants is described in connection 


with the analysis of the cis-compound. 
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The LAME computer program (LAOCOON with Magnetic Equivalence) 


was supplied by its author, C..W. Haigh of the University 


of Swansea, and modified by G. W. Stockton for use on 
the IBM 360/67 computer at the University of Alberta. 
It is similar to UEAITR.17> 
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Literature values for the probable coupling constants 
in cis-perfluoro-2-butene were less useful than they had 
been for the trans-isomer. Two of the coupling constants, 


Ji> and Jo1> were predicted with reasonable confidence 


to be in the ranges -7 to -13 Hz. and +8 to +12 TZ", 


2427-29 The absolute value of J545 was 
24 


but no literature value 


FeSpect ive ly). 
precuctedito tbe Wl Hz, 
for the sign was available. Values for Je_cuc-p (cis) 
in trifluorovinyl compounds range from 24 to 65 Hz., 
with most in the range 30-40 Hee However, the table 
of additive substituent effects “© which predicted the 
value of trans-Jj4 Cogwe Chin Oo 7 eo oy ees) ve eee OG 
cis-J ihe s1¢n ‘of Jr_cuc-p (cis) 1S positive, in most 


— "11° 
compounds, “° but negative in cis-1,2-difluoroethylene-°?>+ 


An attempt was made to visually match the observed 
spectrum using trial-and-error calculations with the 
NUMAR computer program. The magnitude and sign of Ji 
and Jo4 were fairly well known from the ranges given 
above, and the absolute value of their sum N was known 
(by inspection) to be 1.5 Hz. The magnitude of J55 was 
also well approximated; its sign was arbitrarily assumed 


to be positive.! However, the wide range of possible 


t The signs of Jui relative ‘to J55 and of Ji relative 


to Joy affect the appearance of the spectrum, but the 
5 


; 2 
Signs of Jui and J55 relative to Jo and Jo4 do not. 
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Jui values suggested that many trial calculations might 
be needed. It soon became apparent that most of the 
strong lines in the calculated spectra resulted from 

the accidental coincidence of unrelated transitions. 
Small changes in the coupling constants produced large 
changes in the appearance of the spectra, as various 
combinations of lines came into coincidence or separated 
throughout each multiplet. The trial-and-error approach 


was abandoned. 


The only lines in the observed spectrum that could 
be assigned. were those belonging to the doublets of 
spacing N = 1.5 Hz. The intense doublet in the CF, 


multaplet was known to consist of the superposition of 


many lines from various matrices but the corresponding 


We 
doubiet in. the =C-F multiplet seemed likely -to contain 
only the lines originating in the symmetric matrix, with 


all spins F = 4%. 


The low-field (frequency +0.75 Hz. from the =C-F 
chemical shift) member of this doublet was chosen as the 
Mraneaction 1) (see §3.3).to be irradiated in a double 
resonance experiment. It can be seen in Figs. 3-3 and 
5-2 that there are four species X transitions connected 
to one or other of the two coincident species A transitions 
being irradiated. Two of these transitions are members 
Be the species X doublet with spacing N, and the others 


are lines from the Ma = 0 subspectrum. Assuming that 
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FIG, 5-2, INDOR SPECTRA OF THE PERFLUORO-2-BUTENES, 
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the irradiated transitions T are those with frequency 
+4N, the frequencies of the two M,=0 transitions are! 


+4K 2% (K? +12) % (5.2) 


The INDOR double resonance technique!°was chosen, 
because of the high probability that the connected transitions 
would be members of unresolved multiplets. INDOR spectra 
with a poor signal to noise ratio were obtained using 
a modified Varian HA-60I spectrometer (see nye Bey iad Be 
Comparison of a number of repeat spectra indicated that 
the two Ma = 0. transitions connected to T had frequencies 
of approximately +3.4 and -22.2 Hz., relative to the 


-CF, chemical shitt. 


These values were substituted in (5.2) to obtain 
Pies solution K = £18.8 Hz., L = +17.4 Hz. From L and 


Tee, Onan ee OR Seen 


N, Joy and Ji> were calculated to be 
Benew series of trial spectra was then calculated, varying 
the one undetermined parameter, M = J447522> With Ja =a RA, 
Hz. and J55 set sontz che; Calculatedasnectrum closely 
‘resembled the observed spectrum. Other values of Jay 

from -25 to +50 Hz. were used. with appropriate values 


Gi J none of these combinations gave a calculated 


22° 


One cannot determine from the spectrum which of these 


values corresponds to Jj, and which to J,,. 
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spectrum similar to the observed spectrum. 


The INDOR results also provided the remaining 
relative signs of the coupling constants. Of the two 
INDOR lines given by (5.2), the one for which K and Grab ae 
have the same sign is farther from the CF, chemical 
eniitt. ~if this sign were the same as the sign-of N, 
the line would be on the low field side of the chemical 
eee dere 22. 2 Hz). Since: 1t.1s10und, atuo22, 2.2. 
from the CF. chemical shift, when observing at +0.78 Hz. 
meomothe .=-C-F. chemical snulay 2% and N have opposite signs. 


Pees Pens Piven above for the four; coupling constants 


are consistent with this result. 


The INDOR experiment was later repeated using a 
PrukeraiFx-90 "spectrometer. .The resulting spectrum ‘is 
shown in Fig. 5-2. The INDOR frequencies were individually 
Measured under steady state conditions, and had the values 
eer -0.95, +0.75, +3.47 Hz.- from the -CF. chemical 
shift. These frequencies agree quite well with those 
measured in the single-resonance spectra, and confirm 
the original INDOR results. Curiously, the positive 
peak at -0.85 Hz., which was used to optimize the INDOR 


conditions, shows the largest discrepancy (about 0.08 Hz.) 


from the results of the single-resonance spectra. 


The approximate parameters above, with the upper 
set of signs, were used as initial values for an iterative 


calculation with the LAME computer program. The RMS 
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deviation of observed and calculated transitions was 
less than 0.09 Hz. for the parameters given in Table 5-2. 


Observed spectra, and spectra calculated from these parameters, 


are shown in Fig. 5-3. 


An INDOR experiment was also performed (using 
the Bruker“HFX-90). for the trans-isomer. As with the 
cis-isomer, the low field member of the N doublet in 
the =C-F region was observed while a second radiofrequency 
was swept through the -CF. region of the spectrum. The 
Gutecrmost peak, expected to be 140.23 Hz. from the chemical 
shift, was too weak to observe. The three central peaks, 
shown in Fig. 5-2, were as predicted by the previous 
analysis of the single-resonance spectrum. The inner 
peak of the pair given by (5.2) is to low field of the 
chemical shift; by the argument used above, this confirmed 


that K and N have opposite signs for the trans-isomer. 
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aes Discussion 


It was hoped that it would be possible to obtain 
tne absolutersiens, as welJ.as. the magnitudes, of all 
of the coupling constants? in cis-perfluoro-2-butene. 
This has not been possible, because there is insufficient 
evidence to allow definite assignment of the two =C-F 


to -CF, CouuLang ;cCOnNStants +8. 07 and +9.6 Hz. 


The magnitude of each of these values is within 


Picpeanea expected for either the three bond, J. an 
F-C(=)-CF, 
Ocechevrour bond J - 24,27,29 TE seknown that 
F-C=C-CF 
the three bond coupling constant Ji> is negative, and 
28 


the four bond coupling. constant Ja1 IS ,pOsitives big 


Jo oie OFZ. 5, DOK Jui and J55 aie DOS IsGL Me nari Ji> = ~§8.0 


ios, OO th Juq and J55 are negative. 


Evidence from some 2-substituted heptafluoro-2- 


Parenesc 2c! 


tsvcontlicting. int edch cis s1senar: JE-C(=)-CF, 
is smaller in magnitude than JE-C=C-CF," Un slit SaaS Sr, 
Ji> = -8.0 Hz. However, in these same studies, Ju_¢(=)-CF 


in each cis-isomer is at least as large (in magnitude) 


5 


as the corresponding coupling constant in the trans-isomer. 


Since in trans-perfluoro-2-butene this coupling constant 


i -9.1/Hz., one is led to the assignment Ji2 = -9.6 Hz. 


These results can be reconciled in any of three 


ways: 


i) The coupling constants in the perfluoro-2- 
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butenes follow a different pattern from those in the 
heptafluoro-2-butenes,74>2/ 


11) The cis-and trans-isomers were Incorrectly 


identified either in this study or in the previous Se euene | 


iii) The two CF. resonances were incorrectly assigned 


in either the cis-or trans-isomer, in both of the previous 


Mumieee ce! 


The first possibility is quite likely the correct one, 

but is unimportant because it provides no useful information. 
The second can be ruled out immediately by the excellent 
agreement between J55 and Joy values for the trans-~isomers 

in this study and the previous ones. The third possibility 
Belianas caretul scrutiny,. because, if it 18- correct, 


it leads to a clear assignment Ji =t=9),08HZ.. 


here 1S a variety of evidence .. to show that 
in the trans-isomers YF-C=C-CF, Pounediec Jaki oe Thus 
: , : fa | 
any misassignment must be in the cis-isomer. Andreades 
bases his assignment on chemical shift evidence which 
he does not give. Tiers 2 finds chemical shifts of -64.8 


giue-o8.4) ppm. from-CFC1,. for C(1)F, and C(4)F, respectively 


5 
fires -2-chloroheptasluoro-2-butene. ~ Hertherefore 
assumes that the chemical shifts -62.2 and -66.3 ppm. 
in the cis=isomer belong to C(1)F, and C(4)F. respectively. 


This assignment appears to be entirely reasonable. 
The small amount of literature evidence available 


tends to support positive signs for Jaq and Jo: Additive 


substituent epaects 7° predict Jui =k Sc MA eee Oe Clie 


LL. 


git ar esods (wort wrsatseq $e 
¥ 

vitoorTroont 16" 

29s10 stioiveTq oft 4 

(EJQoT7 s9nf stow seo neigest 

~~ onde Teetee 


<8 TOT TOD she yletEL TAPE sas a 
seunsed: Siad tog 


E 
i 
+ 


Cjtecu om 20 Rvorg 


frigaojxOd oy ey «pert vyise7s (bomnts Jud hale ods 


Z 
2" 47kTO rests C2) = 3 ie ect eomnlAyv aP Ie i a 
Cds ‘ yy 
+ ' 5 Sy rt e 7; We to FO FY 
» h } I 4 . H] on 
= : | ; 
P it wor { INIVS 
+” a 
Sirit — ae 
7 = J4 
7 235i 1 bitA 1TOlI0 = ies 13. 40.8 
54 ‘ - i J ee Hey: ’ 
fe 
pa- to 2itinde Legation ce oa . 
! , BY i 
4 F outs 
~~ res” 4 7 ff ° » a 
ytavidosaéex .1(h)0 bos .tGe Las rey 0 work mae 
nwvteayvant. oH os ry Fd -6 H raw yt oa10 oi 
My TeTany oh SATS ss TOU e 2 r 3-8 
i = 


) ; a 
omgy. 2. 90- Daaw: “aan oieriend Lesimeda 03 0 
pt _ eat s 
“is isa5qe2 ot gTEh}9 bite gi(t9 oy guoted 4 ITO & 


owbdsnogesrt vor ttre ‘od GF ae i: 


> 


stdstivve o2aablye, eruse1o? it in, 106 it 


% 


“QV ESTBBA (gol ‘bare’ pk 10? Rit te oui 
ails ot sth 6 yl 


dgie2 


extent that J55 includes a through-space contribution, 


tiamtiay be expected to be positive. >4 


The weight of evidence seems to support the assignment 
given in Jable 5-2. However, further studies to clarify 
the situation are indicated. It would be particularly 
interesting to study cis-perfluoro-2-butene in a liquid 
crystal medium, in order to determine the principle 
components of the J tensors, as well as the signs of the 
isotropic Wed lucs a4 ihe, spectrum would be wieryecompilacated.,, 
but with the aid of the parameters obtained here for 


an asotropic medium, a full. analysis might. be possible. 
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ae The 2,5-Dimethyl-3-Hexenes 


Full analysis of the spectra of cis-and trans- 
2,5-dimethyl-3-hexene was undertaken in cooperation with 
L. W. Kaslander and Dr. F.H.A. Rummens of the University 
of Saskatchewan, Regina, Saskatchewan. The chemical 
shifts and coupling constants, particularly of the methine 
protons, were of interest as.part*of a. continuing study 


of conformations and NMR parameters of olefins.°° 


Cis-and trans-2,5,-dimethyl-3-hexene (Aldrich 
Chemical Co.) were degassed and sealed into 5 mm. O.D. 
thin wall NMR tubes, with 5% tetramethylsilane (TMS) 
added as an internal reference. The infrared spectrum 
(Beckman IR-8) of the cis-isomer had a band at 960 Gta”, 
indicating that a small amount of the trans-isomer was 
Weeseicy, 10 eimpurities were noted ii *the “infrared "Spectrum 
of the trans-isomer. The methyl region of the NMR spectrum 
of each isomer showed weak peaks corresponding to the 
intense methyl doublet of the other isomer. The amount 


of the other isomer present, about 2% in each case, was 


email cenough not to-interfere with the analysis. 


Proton NMR spectra were obtained at 100.0 MHz. 
(23.487 kgauss) with a Varian HA-100-15 spectrometer 
in frequency sweep mode, locked to internal tetramethylsilane 
(TMS). A frequency synthesizer was used to obtain a 
more stable lock frequency (see Appendix E). The spectra 


Were Calibrated as described in Appendix F, yielding 
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line positions estimated to be accurate to +0.02 Hz. 
in the methine regions. It was noted that, for both 
isomers, when a peak was scanned first from the right 


and then from the left. the two traces coincided. 


The NMR spectrum of each isomer can be divided 
into three regions, approximately 5, 2.5, and 1 ppm. 
to low field of TMS, corresponding to olefinic, methine, 


andpme thy laprotons :respectivelys. As a first.approximation, 


each kind of proton was treated as a separate species, 
and the two methyl groups of each isopropyl group were 
assumed to be a single magnetically equivalent group of 
six nuclei. With these assumptions, each isomer could 


be described as an [AKX_ ], Spin system. 


In the discussion below, the olefinic, methine, 
arnetie groupsmareelabelled 1,.-2, and 3 respectively. 
The coupling constant between the two olefinic protons 


is thus J There are two coupling constants connecting, 


ii 
for example, the olefinic and methine protons. These 

are labelled Ji. and Joy? with Ji (indices in increasing 
order) indicating the same-side (close range) coupling, 


and Joy the long-range coupling. 


The general features of the spectra are similar 
for the two isomers. Each methyl region consists of 
an intense doublet with separation about 6.6 Hz., linewidth 
about 1 Hz., and no significant fine structure. Each 
olefinic region has four strong central lines and two 


weak outer lines, all with about 0.7 Hz. linewidth. Aside 
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from a small intensity bias to high field, each olefinic 
spectrum is nearly symmetrical about its chemical shift. 


The methine regions are very complicated. 


These features indicate that the methyl and olefinic 
protons both couple strongly to the methine protons, 
but weakly to each other. Each methyl (region 3) spectrum 
Galawenetueated .as.ithe, X.region. of ja KXe spin system, 
with a doublet splitting equal to Sx: The, 70d eGinac 
(region 1) spectrum is similar to the A region of an 
[AK], spin system with Sexe sma lileeweln: each case athe 
fine structure resulting from small coupling constants 
is unresolved. 

When regions 1 and 3 were analyzed on this basis, 

Hi 


using the [AX] , expressions of PSB’, the approximate 


parameters obtained were: 


Cis trans 
Jui 10.87 15.41 
Ji>> Joy 9.49, -1.04 One ws S 
Ji3 66:62 6.68 


These values were used as the initial parameters for 
trial-and-error calculations on region 2 (methine) with 


the NUMAR computer program. 


The overlapping multiplets of the observed methine 
spectra created some difficulties in comparing observed 


and calculated spectra. Preliminary calculations showed 
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over two thousand transitions in the methine region. 
The peak positions in the observed spectrum had been 
measured very accurately but could not be assumed to 
be the same as the transition frequencies because of 


the influence of other nearby transitions. 


For this reason, Calcomp plots that simulated the 


experimental lineshape were used throughout the analysis. 

It was found that a Lorentzian lineshape, with linewidth 

0.14 Hz. for the cis-isomer and 0.32 Hz. for the trans- 
isomer, provided the best match between observed and calculated 
spectra. The plotter was found to be accurate enough 

that measurement of peak positions in the plots to +0.01 cm. 


was justified. 


Spectra calculated using the [AKX, ], approximation 
gave a poor fit for the central part of each methine 
spectrum, but a reasonably good fit for the outer parts. 
The sign and approximate magnitude of J55 were determined, 
FOrevoOtn isomers, at this stage. [ABX, ] , spectra were 
calculated, but were not used because they were not significantly 
different from [AKX; ], spectra calculated using the same 


parame LErs: 


The final parameter adjustments were made treating 
the spin system as [AXY_], (i.e. olefin protons as one 
species, and methine and methyl together a second species). 
The previously determined parameters were changed very 


little. Peak positions in the methine regions, calculated 
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from the final parameters in Table 5-3, nearly all agree 
with the observed peaks to within experimental error 
(20-02 Hz.); and no deviratron vs greater than 0-1 Hz. 
Figures 5-4 and 5-5 show the observed methine spectra 


for the two isomers, and spectra calculated in the [AXY,], 


Lal ay 


and [AKX, ], approximations using the parameters of Table 5-3. 


Calculated spectra of the methyl and olefin regions (not 


shown) also agree well with the observed spectra. 


The appearance of each spectral region is insensitive 
to small changes in the relative chemical shifts. Thus 
when the coupling constants had been determined, the 
chemical shift for each region was obtained very easily 
from the constant difference between observed and calculated 
transition frequencies. The methine chemical shifts 
BYomestimated ito be accurate’ to +0.02 Hz. (0.0002 ppm.). 
The broader observed lines in the olefinic and methyl 
regions did not allow such accurate determinations; the 
chemical shift errors are estimated to be +0.05 Hz. (0.005 


ppm.) and +0.1 Hz. (0.0001 ppm.) respectively. 


The spectra of the methine regions were surprisingly 
sensitive to both sign and magnitude of all of the coupling 
constants, including Ji3: It is unlikely that the magnitude 
of any of the coupling constants is in error by more 
than +0.05 Hz., and the excellent fit of the observed 
and calculated spectra suggests that most are accurate 


to +0.03 Hz. The small coupling constants (|J|<0.1 Hz.) 
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NMR PARAMETERS OF THE 2,5-DIMETHYL-3-HEXENES 


cis trans 
Chemical Shifts 
(ppm. from TMS) 
=CH (3h) 57, 02 5 713035 
=C-CH (2) a oa 2.2000 
- (CHz), (3) 0.932 0.950 
Coupling Constants| 
OH Zs) 
Jui 10.88 15.40 
Jyo> Joy OPA Tee LZ OS Oy a0 
J55 0.39 0.60 
Ji3: Jay O, 208 0.8 -1.08 ,f 0.0 
Jos Jz, 6.59 51°0.0 el Pare Wise OE 


t In pairs such as Ji and J51> Joy is the long range 


coupling constant. 
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FIG, 5-4, OBSERVED AND CALCULATED SPECTRA, 


C=C 
\ 
H —-HC(CHs), 


whl fl ! 
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FIG, 5-5, OBSERVED AND CALCULATED SPECTRA. 
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were varied in smaller steps; their values are probably 
accurate (0 40.02 Hz. Thecsix bond coupling constant 
Jz5 in the trans-isomer was varied from -0.1 to +0.1; 
the best agreement between observed and calculated 
intensities was obtained for Jz =et0.028 BThvs wake 


1s so near zero that it may merely be compensating for 


bio giercordereeffects ignored in. the [AXY J, cabculation. 


Piet oeObVlLOUSs that, tor dil practicdaye purposes), 
the [AXY_], approximation adequately describes both of 
these spin systems. The few residual differences between 
Observed and calculated spectra may equally likely be 
the result of small errors in the coupling constants 
or of the assumption that the olefin protons can be treated 
aseauseparate species. There is no- evidence torsuggest 
that the six methyl protons of each isopropyl group are 


not magnetically equivalent. 


One puzzling minor observation should be noted. 
The linewidth needed to adequately match the spectrum 
of the trans-isomer is more than twice that needed for 
the cis-isomer. Higher-temperature spectra, which have 
not yet been fully analyzed, show noticeably sharper 
lines for the trans-isomer above 458Gea efheresastnor 
sufficient evidence to indicate with any certainty the 
cause of the broadening, but it could perhaps be the 
result of high viscosity in the neat liquid at ambient 


temperatures. 
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Dre. Trifluoromethyl Phosphorus Compounds 


The NMR spectra of a series of symmetric compounds 
containing -P(CF,), were analyzed as [AX], spin systems, 
in cooperation with Drs, Ass Pinkeriton andt RiméLn Cavell 
Ofte this Univensi ty<ie The analysis for (FC) ,P-NH-P-(CF2)., 
a typical member of the series, is presented briefly 
Dorceasea hurther illustration ofthe ease with which 
the spectra of large, highly symmetrical spin systems 


can be calculated using the methods described in earlier 


Chapters of this Thesis. 


The [ (CF) ,P],NH was prepared and purified by 
Peeertikeorco, Using, the usual procedure. >“ A sample 


was dissolved in CCl, with® added» CFC1. as an internal 


3 
reference and sealed in a 5 mm. O.D. NMR sample tube. 
19F spectra were obtained using a Varian HA-100-12 
spectrometer operating at 94.1 MHz. (23.487 kgauss). 

A frequency synthesizer was used to allow large offsets 


from the CFCl, lock (see Appendix E). *'P spectra were 


3 
not obtained. 


Linewidths were somewhat greater than in the other 
compounds studied. This may have been the result of 
coupling either to '*N, or to the NH proton. Since it 
produced no resolvable fine structure, the NH proton was 


ignored in the calculation. 


The phosphorus and fluorine nuclei of this compound 
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= J 


labelled species 1 and species 2, and the three-bond 


four-bond P-F coupling constants are labelled J 


Bub 
Joy respectively. One coupling constant combination, 
127924 = 91.2 Hz., is available by inspection from 


separation of the intense doublet that dominates 


PE spectrum (Fig. 5-6). The ?°F chemical shift 


Pomoogoeiz. to'high field of CFC1, C5 So eed Dili) 


with the assumption J55 = 0, using the procedure of Harris. 


The separation of the innermost line outside the N doublet 


iat uuadanva lucsnot tJ and L = Ji27J were obtained, 


11 21 


from the most intense line inside the N doublet on the 


same side gave Jy -2254.0 iZseelnds "Valucewasetnen 


substituted in Harris' equation [5], which gives the 


£25 


14 


separation of the innermost pair of lines as (Laan) ce ee |e 


Fromethe result, L =.82.2 Hz., and the»known. value of \N, 


Ji 


and Joy were piven as 756.7 and 4.5 Hz. respectively. 


These values were used as starting parameters for 


trial-and-error calculations with the NUMAR computer 


program. The best visual match between observed and 


calculated spectra was obtained for Gineis.s) 


UA 228.0 
ion din 86.9, 4.3 
ing 0.7 


All of these coupling constants have estimated errors 


of +0.1 Hz. The spectrum calculated using these values 
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is shown in Fig. 5-6. Observed Spectra are shown above 
the corresponding calculated spectra with the central 
lines and the outermost lines in an expanded scale. 
Frequency scales below the calculated spectra are in 


Hz. with an arbitrary zero. 


The signs of Jui and J55 are the same, and the 
Signs of Ji and Joy are” the samé;. the signs of Jui and 
Joo relative -to Jy and Joy do not affect the appearance 
25 a : F 
Of the spectrum. The positive signs given are consistent 
with literature eaies. 2 but should be confirmed by a 


double resonance experiment. 


Other compounds whose NMR spectra have been analyzed 


in a Similar way include [(CF.)..P(=S)],0, LOGE Py,0, 


3)2 
[(CFz),PJ.S, and [(CF,),P(=S)S-],. In most cases’ the 
spectra could have been calculated by perturbation necnotnoe 
but exact calculations with NUMAR program were much faster 
(ten seconds computer time; typically, on an IBM (360/67 


computer), and provided plots that allowed a more accurate 


comparison of observed and calculated spectra. 
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CHAPTER 6 . CONCLUSIONS 

A new approach to the calculation of the NMR spectra 
of molecules with twofold frame Symmetry has been developed. 
A total spin operator A is defined for each symmetrically 
equivalent pair of magnetically equivalent groups 

A (pair 2) = B (2) + F ce’), (6.1a) 
and for each invariant magnetically equivalent group 

A (group i) = ‘Sule | (6.1b) 
For each pair, an antisymmetric operator B is also defined: 


beinosret).= © (2) ~ F (2°). (6.2) 


Mhestwoudistinc&® coupling constants J and J ' 
near fat 

connecting the magnetically equivalent groups of pair 2 

to those of pair m, are replaced by the symmetric and 


antisymmetric combinations 
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The nuclear spin Hamiltonian for the molecule is then 


written 
A 2 vil) + PI pp FO) FO) 


+ 2 J4;hG) AG) | (6.4) 


1 <j 


+ Jmeb(£) *B(m) , 
m 


In (6.4), sums over i, j, ... include both invariant 
groups and pairs, while sums over 2, m, ... include only 


the pairs. 


A suitable basis for this Hamiltonian is constructed 
from eigenfunctions of A? and Boy for. allthe anya ie and 


groups and pairs: 


} | A, M, >. Cor) 
When F # F” for any of the symmetrically equivalent pairs, 
these product functions are not eigenfunctions of the 
molecular twofold symmetry operator C,. When F = F” 

for all pairs, the symmetry of each product function 


1S given by 


C, I A, ,M;> = Cig) eee aes) Th SAM (6.6) 
pee tira cont £ i of 
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The basis functions belonging to the symmetric (A) repre- 
sentation of C, form one matrix, , and those belonging 
to’ the antisymmetric (B) representation form another. 
Thus this basis yields all of the symmetry factoring 

that results from twofold frame symmetry as well as 


magnetic equivalence. 


Explicit expressions have been derived for all 
Pee tiesiatirix elements or the Hamiltonian (6°.4)° 17% the 
basis (6.5). Although these expressions can be used 
for hand calculations, their chief advantage 1s their 
suitability for use in a general computer program. Unlike 
previous expressions for molecules with frame symmetry, 
they are’ valid for molecules of arbitrary size, with 
either no frame symmetry, or twofold frame symmetry. 
(For molecules without frame symmetry, (6.4) and (6.5) 
reduce to the corresponding equations of the Composite 


Particle method. % 


This symmetry factoring has been combined with 
‘factoring by the z components of molecular total spin M, 
and species total spins Mo» to provide a systematic, three 
stage approach to the factoring of the Hamiltonian matrix. 
It has been shown that this approach to factoring leads 
in a particularly straightforward fashion to the selection 
rules governing transitions between Hamiltonian matrices, 
and to a description of each species’ spectrum in” terls 


of subspectra connecting separate series of matrices,. 
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NUMAR, a flexible and efficient computer program, 
has been written to allow the use of all three Stages 
of factoring for calculations on large spin systems. 
The user has the option of making quick initial calculations 
assuming a high degree of symmetry and X factoring, followed 
by more accurate calculations eliminating any factoring 
that has proved unjustified. Quantum numbers are printed 
with the energy levels and the transitions to simplify 
iieeiventiti1cation of subspectra.* Two plotting routines 
are included in the program. The one provides bar plots 
on the computer's printer. for immediate use in the early 
Stages of analysis. The other provides a variety of 
lineshape and scaling options to allow accurate comparison 


of the fine details of observed and calculated spectra. 


The new factoring procedure has been applied in 
the analysis of the NMR spectra of a number of molecules 
with twofold frame symmetry. Preliminary parameters 
for the perfluoro-2-butenes were obtained from the literature, 
and from InDoR! © spectra. Preliminary parameters for the 
2,5-dimethyl-3-hexenes were obtained from approximate 
analyses of their methyl and olefin regions. Preliminary 


parameters for [ (F.C) ,P],NH were obtained from Harris!?° 


He In each case, the preliminary parameters 


equations. 
were refined by a series of trial and error calculations 
with the NUMAR computer program; the resulting final 


parameters gave calculated spectra that matched the observed 


spectra within experimental error. 
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The approach to NMR symmetry factoring developed 
in this Thesis is the most general presently available 
for calculations on large spin systems. It suggests 
several opportunities for future research, particularly 
the development of a similar procedure for molecules 
with Coy or D. frame symmetry, and the extension of the 


equations and computer program-to include partially 


oriented molecules in liquid crystal media. 
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APPENDIX A. SYMBOLS AND NOTATION. 


Spin Systems 


foe es Bee Nuclei with similar chemical shifts; 


One ySpeg Les. 


Ao 0 Re Nuclei widely separated in chemical 
| SHaL Cet LOmen bole ano tier 1s peGles. 
An or Xn Magnetically equivalent group of n 
nuclei. 
[.--], Indicates an n-fold element of frame 


SVIMEEnY. ssl nesnotallon ror da ty pricar 
member of each set of nuclei (or 

of magnetically equivalent groups), 
related by the symmetry element, 


is enclosed in the brackets. 


For example, p-fluorotoluene, at 24 kgauss, can be treated 
as an Az[KL],X SUlIeSy SLM owl ClimitikLe ee speGles sen lit. itie 
other commonly used notation, it would be called 


AZKK box. 


Chemical Shifts and Coupling Constants 
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V5 Chemical Shift i, in Hz., relative 
to Vo of the species containing group 
OTe pains, 

Jij Coupling Constant, in Hz., between 
Magnetically equivalent groups or 
pairs. 

Sng? Jap The short range and long range coupling 
constants connecting pairs p and q. 

Jop- The coupling constant between the 
two magnetically equivalent groups 
OLapa binges 

Jom? Jnl Symmetrized coupling constants: 

J 2 4..(J +) 
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Aldeespin operators used an this Thesis are 
dimensionless. 
Operator 
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Spin operator for a single nucleus. 
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AS 


equivalent groups of nuclei. Identical 


to ig for invariant groups. 


Species S total spin operator (z 
component. 


Antisymmetric operator, complementary 
to ike for each symmetrically equivalent 
pair of magnetically equivalent groups. 


Not defined for invariant groups. 
Hamiltonian operator. 


Molecular transition operator. 
Transition intensities are proportional 
to the absolute squares of its matrix 


elements. 


The corresponding operator is given to the right 


of each eigenvalue. 


A (At+1) 
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A, (S) 


Molecular total spin (z component) 
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Basis Functions 


In kets of the form |F,m> the first eigenvalue 
always corresponds to the squared total spin operator 
2 


F , and the second to the corresponding z component 


operator set 


F? |F,m> 


E(PtLy sm 


Be |F ,m> m |F,m> 


The |I,m> and |A,M> kets also follow this convention. 
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APPENDIX B. THE NUMAR PROGRAM (FORTRAN LISTING). 


NUMAR is written in FORTRAN IV at the level of 
the IBM FORTRAN G and FORTRAN H compilers. It has been 
tested on the IBM 360/67 computer at the University of 


Alberta under the 0S/360 and MTS operating systems. 


The routines TIME and CREPLY, used in INPUT and 
BPLOT, are specifics to the MTS operating system. They 
may be omitted or replaced without significant effect 


on the remainder of the program. 


Double-precision (REALs8) arithmetic is used in 
HAMINs “to reduce rounding. errors, and half-precision 
(INTEGER*2) integers are used for some quantum numbers 
to save memory space. Neither of these choices is essential 


to the operation of the program. 


NUMAR is designed for efficient operation as an 
overlay (or chained) program. The small main program, 
beth Common blocks, and the FORTRAN library form the 
base segment. The single overlay region is used in turn 
for each subroutine called by the main program (and any 
subroutines it calls). NUMAR as presently dimensioned 
can execute in a storage region of less than 100K bytes 


(25000 words). 


See also Appendix D, Instructions for NUMAR. 
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APPENDIX C COMPUTER REPRESENTATIONS OF SPIN QUANTUM NUMBERS 


Since spin quantum numbers are multiples of 0.5, 
it seems reasonable to represent them in a computer by 
floating point numbers.' Unfortunately, different sequences 
of floating point arithmetic operations that should yield 
the same result often do not because of rounding errors. 
Logic based on floating point quantum numbers is thus 
unreliable. A suitable internal representation for a 
spin quantum number is an integer equal to twice the 
true value. All equations for matrix elements must be 


modified accordingly. 


Memory space can be saved by 'packing' a set of 
related quantum numbers, such as the A values of all 
groups. Each group's A value is multiplied by the corres- 
ponding element of a vector of multipliers derived from 
the maximum A values for the molecule; the sum is a unique 


iaeevervrepresentingr thes entires sete 


Bach multiplier, starting: from: the right, is the 


1 Non-integral numbers are normally represented in a computer 


by the product of a fractional part and an exponent, 
together known as a floating point number. The principle 
is the same as that of standard scientific notation, but 
base 2 rather than base 10 is normally used. Many numbers 
that are exact in scientific notation become endless 
fractions in floating point, causing unexpected rounding 


errors. 
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product: 

(maximum A value next right, +1) x (multiplier next right) 
For example, the maximum A values of an A[BX¢], Spin 
oretarcogl) 0, or, an internad representation, 1 2 12. 
PUteespilesystem's multiplier vector.is 39°13 1, and the 
Beto. Avvalues 4% 0 5 is packed as’ 49. To unpack, one 
reverses the packing process: 49+39 = 1, remainder 10; 


10:13 = 0, 


Unpacking is much more complicated if negative 
numbers are included in the set being unpacked. If M 
(or Mg) Walucsmeare packed sitheveshould stirs t. be wmade 
non-negative by adding to each a number equal to the 
largest possible absolute value of M (or Mg) POtepi at 
peoup Onrespecies.. .In.an A,B xX system othe =1imi towtor 
Mo Values of the two species are +2.and +3/2; in internal 
Poti and +s, a betore packing +4 -andu+3 srespectively 
are added to the Mg Values! -OLVin ves aneese cd. sure COwss 
Shas erouv0cuw l(husathe packed representation sot Mo values 
Ottmann. 1S 58. 

Packing is particularly advantageous when the 
packed representation of a set of quantum numbers can 
be used directly, without unpacking. As is pointed out 
Sie eo basis functions can be sorted ito matrices, 
by their packed Mo values, and packed A and Mg values 
can be used directly to determine transition selection 


rules. Since there is no direct use for packed group M 


values, they are best left unpacked. 
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APPENDIX D. INSTRUCTIONS FOR NUMART 


. NUMAR is a new program for the calculation of 
NMR spectra from chemical shifts and coupling constants. 
It makes full use of the factoring that results from 
magnetic equivalence, twofold symmetry, and large chemical 


Snitt ditterences. 


SOME DEFINITIONS 


ool bo enue TOPIC SPECIES: 

Usually this is just what it seems. Protons in 
a molecule are one species, ‘°F nuclei another, and so 
PoceNUCiScinoflithée samerisotopacispeciesebut-farsapart 
in chemical shift may also be entered as different species. 
PrommthisevxX-approximation"” to be valid, all coupling 
constants between nucle1 in different species must be 


much smaller than the chemical shift separations. 


CASE: 
The set of data cards necessary for calculations 
on a molecule is a CASE. Each computer run may include 


calculations on any number of cases. 


: University of Alberta Nuclear Magnetic Resonance. 
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Here 
The cards making up a CASE may be grouped into 
several TYPES such as chemical shift cards, coupling 


constant cards, etc. Each TYPE has its own layout of 


columns for punching. 


GROUP: 
In these instructions a GROUP always refers to 


the collection of nuclei described by a single TYPE 4 


card. A GROUP may consist of: 


a single nucleus of spin % up to 9/2. 

- a magnetically equivalent group of up to 9 nuclei 
OFnSDiINnAweE 

- a chemically equivalent pair of either of the 

above. 


The groups*are numbered (for coupling. constants etc.) 

Seauentially from the first group of. the first isotopic 
Species to the last group of the last isotopic species. 
A chemically equivalent pair entered on one card counts 


as one group not two. 


(There is just one time that you count two for 
a chemically equivalent pair. When adding up groups to 
see if the molecule will fit within the program limit 
of nine groups, each pair counts as two groups. This 
is the result of an internal requirement in the program 


and does not affect the numbering of the groups.) 
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ITEM: 

An ITEM is the piece of data punched in any single 
region of a data card. For example, one data ITEM on 
a TYPE 4 card is the chemical shift, which is punched 


in card columns 10-19 (inclusive). 


ORDER OF TYPES 
A typical case includes up to six different TYPES 
of card: 


1) One TITLE card. 
2) One or more COUPLING CONSTANT cards. 


For each species, a-set consisting of: 


a),Une. LSOFORIC.SPECIES card. 
4) One CHEMICAL SHIFT card for each group of nuclei 
insthis rsotoprc“specres. 
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Nrevettibnke ARE SEVERAL ISOTOPIC SPECIES the order of types 


roi 
Nees tet i) sl ois eel as 


Square brackets are used to show that some types, or 


sequences of types, will not appear in every CASE. 


FORMAT OF DATA CARDS 
In the detailed description of each TYPE only 
the main items on each card are included. For less common 


options see the separate sections at the end. 


A format code lietter and a card column range are 
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given for each data item. Integer (1) values must be 
numbers with no decimal point, and must be punched as 

far right as possible in their column range. Floating 
point (F) numbers must have a decimal point. For both 

I and F, a blank and a punched zero are treated identically. 


Alphabetic (A) data may include any normal keypunch character. 
Sometimes examples are included in parentheses. 


In laying out the cards an attempt was made to 
have similar kinds of data appear in similar placed for 
eecn 1 YPE. 

~ecne- letter Loin column lo ofvascardeaiways=ana1cates 
the last member of a group of cards of the same type. 

pecvetyeinteger "(Ll )rnumbereend seein column e720 tes 

- every floating point (F) number occupies either 
10 columns beginning in column 10, 20, etc. or 5 columns 
ieee ive. COLUMN 205. °25.,.-50, 1.55.4, 6bC. 


=-.every title begins in column 8. 


DEFAULT VALUES 

A default value is assigned to many data items 
if nothing is punched in the card columns. For example, 
if column 2 of a TYPE 3 card is left blank, it is assumed 


that the isotopic species has spin %. 


When the item is one of a pair of frequencies 
that define a frequency range, the default value is assigned 
only if both of the frequencies that define the range 


ane equal. This will be the case if both are left blank 
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or 0.0 is punched for both. Equal positive or negative 
frequencies have a special meaning for some frequency 
ranges, as described on the following pages. Note that 
a~frequency range of blank to 100.0 is treated as 0-0 
to 100.0 and does not give a default value for the left 


hand limit. 


NORMAL INPUT FOR NMR PROGRAM NUMAR 


The most important, items, are indicated by + in the margin. 
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> Al Must always be an asterisk *. 
Ana? Normally blank. Nonblank only if this 
case is one of a series with Short Form 
Input (see p.D-9) 


AP 3 Normally blank. Master Output Control 
(see p.D12). 


Ae 8-79 Anything that you want for a title. 


19 MAS oafde COUPLING CONSTANTS 


oe Al Themietterebeif this as they last. coupling 
constant card for this case; otherwise blank. 

=e Led The first group number. 

a I 4 The second group number (normally larger). 


(lo 2)enrore neu orhexr i transgs Coupling 
between two chemically equivalent pairs, 
reverse the order (2 1). Repeat the number 
(2 2) for the coupling between the two 
members of a pair. 
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TYPE s 


- F 10-19 


ly EWA) 


Pie 54 


F 40-49 


| Ee 


Continued 


The coupling constant in Hz. 


ISOTOPIC SPECIES 


The letter(Lyif this*:isethe lasty isotopic 
Species for this case; otherwise blank. 


An integer equal to twice the nuclear 
spin of this species. The default (if 
left blank) is 1, meaning spin %. 


Normally blank. Species Transition Table 
Control (see p. D-12). 


Blank if no plots are to be made for this 
species. If any nonblank character is 
punched}$*one or more TYPENS cardspi{one 
per plot region) must follow the TYPE 4 
cards for this species. 


A name for this species. 


Lower intensity limit for keeping transitions 


calculated for this species. Default 
(if left blank) is 1/10000 of species 
LOCal intensity. “A nevative value cancels 


all Craisdt Lon, Calculations  foretni ls (species, 


Spectrometer frequency in MHz. (Used 
only; to calculate ppm values for printing). 
Default 60.0 MHz. 


Grouping width. wit a. Condensed frequency 
table has been selected for this species, 
this determines 1ts,resolvtion, 7 in ail 


transition tablés a blank liné 1s inserted 
at intervals of.100 times the grouping 
Witte wee oUt ie U Om iaee 


Upper frequency limit for this species' 
transitions. 


Similar lower limit. Frequencies outside 
the range are not stored at all, so are 
missing from plots as well as tables. 
Default (if columns 30-49 are blank): 

no frequency timitssfortthesspeciestsedo 
not confuse this option with the plot 
frequency range on a TYPE 5 card, which 
affects ‘only that. pilot. 
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F 40-44 


CHEMICAL SHIFTS 


The letter L “Ei this isthe las t*chemical 
Sit te-to0 retiis Species; jotherwise billank- 


Normally blank. If this card is for a 
chemically equivalent pair of groups, 
GHewGianacters 2. 


The number of nuclei in the magnetically 
equivalent group. (A chemically equivalent 
pair of methyl groups would be entered 
as*2*3 in columns 2-4). 


The chemical shift in Hz. from the species 
reference. 


PLOTS 


The*letter “Lei fethis»1rsetheslaseewiot 
card fonathis species, yotherwisesblank: 


Plot tyne a Blank sor zerouclves-aglaine 
printer bar plot... Other values give ‘Calcomp 
plots with different lineshape” functions: 

ye Gauss an sec bOVeNn C24 ane 5.4 por 4h 
modified Lorentzian with increasing fourth 
power contribution. . for most» Calcomp 
PEOESMUSE 7). 


Urdinarily ignored iby program. For plot 
cards-in-Shortsrorm Input tp. p-9)),) the 
species. 


Wo title-tor sthis plot regiom. 
The left. frequency limit for*this plot» region. 


The, right limit. 1fsthe two slimes sare 
equal, all transitions of *thiSespecies 

are included. Equal positive limits give 

a plot in increasing frequency order; 

equal negative (or zero) limits give 
decreasing order. » In the bansnlot, sany 
large regions of blank baseline are skipped. 


Horizontal Scale. For. bar plots, the 
width (Hz.)}, ofeach bar (default 035). 
For Calcomp plots: if positive, Hz/cm.; 
if negative, the absolute value is the 
plot width in cm.; default 20 cm. plot 
width. 
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TYBEMS Continued 


F 45-49 Vertical Scale. For all plots, the default 
is full scale for the most intense line, 
and any positive value is used to multiply 
this scale. If the scale is negative, 
its absolute value is the scale factor 
in symbols (bar plot) or cm. (Calcomp 
plot) per intensity unit. 


PeSU- oA Left partial width at half height (Calcomp 
plot only). 


LS re) Rightupartiaiswictho. [herre summers tne 
Line shahiwtd th. wi(default Oaks sbs0 fhowe 0.30). 


SAMPLE DATA 


iiertirst, thirty columus 0: eae ctypicais set Ul-data 
for 1,l1-difluoroethylene are shown. In this example 
the remaining columns (31-80) of all the cards would be 
Drank. To the right.of the. vertical lines are card: TYPE 
numbers and comments, none of which would appear on actual 


data ‘cards: 


- 1,1-DIFLUOROETHYLENE Bi ee ad eg od = 
J 
1 4.8 Z Jon 
Zz G.7 2 Jur cis 
1 a Picges. 2 Jur trans 
L22 36.4 2 Jer 
# PROTON 3 isotopic species 1 
(no plot, condensed 
table) 
[is opal ru 4 Chemical Shift of 
group l 
L P FLUORINE SeesSOLOPIC. SNE CLES «2 
(plots) 
i2 0.0 4 Chemical Shift of 
group 2 
L COMPLETE Se Plots cegion w 


(using defaults) 
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SHORT FORM INPUT 


Frequently several cases in a run are calculations 
on the same molecule with changes in chemical shifts, 
coupling constants, etc. The first case of each such 
series 1s exactly the same as it would be for a single 
separate calculation. For the others, a special short 
form of input is available in which only changes from 


the immediately preceding case are given. 


Each case for which short form input is used 
DewiiiseWith a 1YPE-1 card that nas a nonblank*character 
in column 2. This is followed by one or more modified 
TYPE 2 cards. There are three kinds, labelled 2A, 2B, 
andecG Delow; they may be intermixed in,any order. Only 
Picmlascwlrrs 2 Card, WNethier 10-15 2A su2by.0r cGy nas 
Merrecolunniel.. [here are no: TYPhsjeor ir bede cards, 


Bue there Can be moditied TYPE’ 5 cards. 


PYee 2A. COUPLING CONSTANTS. Identical to a normal 


Bype 2 card. 


TYPH.ZB. .~+CHEMLCALSHIFTS..,Column.Z;; .the,groupsorypair 
number, in the same numbering scheme used for coupling 
constants. Column 4 must be blank. Columns 10-19: 


the, chemical shift. 


TYPE 2C. SPECIES DATA. Column 2 must be blank. Column 4: 
the species number. Columns 3 and 20-59: the same 


as the corresponding columns of a TYPE 3 card. 


In each case with short form input, any coupling 
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constant or chemical shift for which no new value is 
given remains as it was in the immediately preceding 
case. The default for an item left blank on a type 2C 
card is whatever value was in effect in the immediately 


preceding case, whether or not that value was itself 


ace tault, 


Ordinarily all plot regions previously specified 
tomalieineet rect ‘fora case with short, form input. 
Particularly when default frequency ranges have been 
used on.the type S cards this is: likely to be satisfactory. 
To set new plot ranges, use the letter P as the nonblank 
character in column 2 of the TYPE.1 card. ALL previous 
plot nanges ane then cancelled, and at Least one TYPE 5 
cand must follow the type 2 cards of the case. The number 
of the species to which the region belongs must be placed 
in column 4 of each of these type 5 cards. (If column 4 
is left blank, the region belongs to the nonexistent 
species 0. A single species 0 card can be supplied when 
the intention is to cancel all plots for this case.) 

The new plot regions replace the previous ones as the 
default plot regions for following cases with short form 
input. 

The sample short form data cards below could be 
used immediately after the cards on p. D-8 to cause two 
more calculations on 1,1-difluoroethylene. Fluorine 


transition calculations are cancelled, and the output 
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is reduced to a condensed proton transition table (see 


Controlling Output). 


A ei 1,1-DIFLUOROETHYLENE, CONTINUED. 
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CONTROLLING OUTPUT 


The printed output of NUMAR consists of a series 
Giabies,. most of which. can be omitted if desired. - The 


desoult: 3S.to print all tables: 


INPUT DATA. Group sizes, chemical shifts, coupling constants, 


etc. [his table is always: printed: 


SPECIES TRANSITION LIMITS. These few lines that follow 
the input data are not printed for cases with short form 
input. 

ENERGY LEVELS. Includes quantum numbers useful in tracing 
subspectra. Optional. 

KEY. Gives all transitions that appear in the following 


transition tables, arranged by transition number in groups. 


Shows which energy levels a transition connects. Optional. 


TRANSITIONS. Each species has its own optional transition 
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table. A full transition table includes the line number, 
frequency, intensity, and subspectral quantum numbers 

of all transitions of the species. In a condensed transition 
table, each series of transitions whose successive frequencies 
are separated by less than the Grouping Width (TYPE 3 

card) is combined; only the intensity-weighted mean frequencies 


and the corresponding total intensities are printed. 


The Master-Output Control (TYPE<l.card} indicates 
which of the optional tables are to be printéeder-LGijit 
is blank, all tables are printed, including full transition 
tables for all species. Some nonblank characters, listed 
below, have special meanings; all others cancel the Energy 
Level and Key tables, and give condensed transition tables. 


Thesspeci ai, control characters are: 


E: Energy Level table, no Key; condensed transition 


tables. 


K: Key, no Energy Level table; condensed transition 


tables. 


F: Full transition tables; no Energy Level table or 
Key. 
X: No Energy Level, Key, or transition tables. Used 


when only plots are of interest. 


If the Species Transition Table Control (TYPE 3 
card) is left blank, that species' transition table output 
is controlled by the Master Output Control as described 


above. If any nonblank character is punched it overrides 
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the Master Output Control for that species only. A Species 
Transition Table Control has no effect on the Energy 
Level or Key tables, but otherwise different characters 


have the same meaning as above. 


These Output Controls affect tables, but do not 


affect plots. 


In cases with Short Form Input, the new Master 
Output Control on the TYPE 1 card overrides any Species 
Transition Table Control characters from the preceding 
cases. To override the Master Output Control, include 
a TYPE 2C card with a suitable nonblank character in 


column 3. 


One can completely eliminate a species, or restrict 
its Output to. a. limited frequency range; See? columns 20-25 
aieero-59 Of the TYPE) 3 card. , Ihésesoptions ariecu.plots 
aseweill aS tablesw If a species" transitions have been 
completely eliminated, any plot regions specified for 


it are ignored by the program. 
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HINTS 


Since there is a separate transition table for 
each species, one has considerable freedom in entering 
chemical shifts. For example, when all the nuclei of 
a species have the same chemical shift, symmetries in 
the transition pattern can be made more obvious by using 


Uru ase tne: shift. 


Transition frequency calculations are made to 
seven significant figures. Ordinarily this lets you 
use the actual shifts from the species reference with 
accuracy of better than 0.01 Hz. However when all the 
shifts of a species are several kHz. to one side of the 
reference (as in many fluorine spectra) it is a good idea 
to add a- constant to all the shifts of the species, reducing 
their absolute values. This amounts to a change of reference 


frequency. 


Be especially careful about plot cards. There 
must be a nonblank symbol punched in column 5 of the 
type secard 1f there are any plot (1TYPEo5) cards for 
fhetespecies. Conversely, if-column9> 0f a TYPE 3) card 


is blank there must not be any plot cards for that species. 


Some current limits: 9 groups, 9 plot cards, 
6 species. There are several internal limits that depend 
on the number and size of groups. If you exceed one, 
greater use of the X-approximation (species factoring) 


may help. Maximum use of available factoring also saves 
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considerable computer time. 


Occasionally it is necessary to use the species 
frequency limits, and perhaps raise the minimum intensity.. 
A very large molecule may excede the limit of 2500 transitions 
perespecics, making 1t impossible to sort them all. “One 
must recalculate the spectrum with a small enough frequency 


range that under 2500 transitions are included. 


The remaining transitions of interest are then 
calculated in the immediately following case or cases, 
using new frequency limits specified on TYPE 2C cards 
(see Short Form Input). The program saves considerable 
computer time by going straight to the transition calculations 
when it finds a case with Short Form Input and no changes 


in chemical shifts or coupling constants. 


If fairly, weak lines are of interest,. one must 
be careful to specify a sufficiently low Minimum Intensity 
Pimit. on the ype 3. card. The default valuesmay give 
a calculated spectrum in which some weak lines are missing, 
and others (which are superpositions of many weak transitions, 
some below the limit) have less than their correct relative 


intensity. 
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APPENDIX E. ITERATIVE ADJUSTMENT OF NMR PARAMETERS 


Several methods of obtaining chemical shifts and 
coupling constants have been discussed above, particularly 
in 53.3 and §4.1. None of. them yields exact values for 
the parameters. This is obvious when, for example, 
borderline X factoring is used deliberately to obtain 
approximate values. It is also true when parameters 
are obtained from the algebraic solution for a simple 
subspectrum, since the experimental frequencies used 
in the calculation are not known exactly. The influence 
of experimental error can be greatly reduced if the observed 
Piequencies of other Lines cineadditioneco whe 1tewe that 
belong to simple subspectra, can be used to help determine 
the values of the NMR parameters. 


Various least squares methods?’ >18 


have been developed 
in an attempt to obtain the best available parameter 
estimates. Each method uses a computer program, and 

invOLVes: a ‘Similar series of steps, The program calculates 
aespectrum from a set of starting parameters previously 
Obtained. Ihe analyst then matchesvup lines in this spectrum 
with their observed counterparts, so that whenever possible 
an observed frequency is assigned to each calculated 
transition. The program then adjusts the approximate 


chemical shifts and coupling constants to minimize the 


sum of the squares of the differences between observed 
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and calculated frequencies. 


As the relationship between the NMR parameters 
and the transition frequencies is not linear, the corrected 
parameters are themselves only approximations CO rie 
true Teast "squares ‘parameters: * "However, if the original 
assignment was correct, they will normally be closer to 
the least squares values than were the starting values. 
Thuswatter several repetitions, of the adjustment, step 
one. obtains a set of parameters that represents the, best 
available estimate of the experimental parameters consistent 
with the assignment made. 


A useful byproduct of the iteration procedure 


17,18 is an estimate 


in all of the common iterative programs 
of the probable error in each NMR parameter, derived 
from the residual differences between the observed and 


Galculateduspectra. 


All of the programs mentioned above lead to similar 
values for the least squares parameter estimates, but 
they make use of two fundamentally different approaches 


Portie cCauctitation. 


In the eS based on the NMREN/NMRIT programs 
of Swalen and RITE © attention is focussed on the 
pattern of energy levels of the spin system. Each observed 
transition that is correctly assigned provides information 
about the relative positions of two energy levels; collectively. 


all of the observed transitions assigned to a matrix, 
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must provide information about the positions of all of 


the energy levels of the matrix In the NMREN step, 


5° 
this information is used to provide a linear least squares 
solution for all of the 'observed' energy levels, E obs? 
with each matrix, treated separately. 

In the NMRIT step, an approximate Hamiltonian 
matrix Ha pprox is calculated (using the same starting 
parameters used previously when assigning the observed 
transitions to the energy level diagram) and diagonalized. 

=a) . 
= ie 
approx i: Eealc We aRe 


PieecransrormatLon matrix U that relates the Dasis and 


diagonal representations of H is used to generate 


approx 
from Bee the diagonal matrix elements Hes of an improved 


Hamiltonian matrix. 


Cian 


Thesiiudiagonal matrix eYenénts! tare-used in adbeast squares 
solution for a new set of NMR parameters. An ingenious 
feature of this second least squares step is that no 
explicit matrix inversion is required, because the normal 

: ‘ 18a : d 
equations are diagonal. The new parameters are use 

; MRIT is repeated 

to calculate a new Ha pprox? and so on N 1 p 


until additional iterations fail to reduce the residual 


errors in the Has values. 
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The NMREN/NMRIT approach is not well suited to 
calculations with full factoring. In the NMREN step, 
chains of transitions extending throughout a matrix. 
are needed to establish the relative positions of the 
energy levels. As a result, the least squares solution 


FOrne of each matrix, requires the inversion of a 


obs 
large coefficient matrix equal in size to the matrix. 
Even if sufficient main storage is available for NMREN, 
other difficulties remain. Only the diagonal matrix 
elements are used in NMRIT to calculate new parameter 
estimates. If the frame symmetry basis (2.24) is used, 
the symmetrized coupling constants Jnl of aba am(2 922) 


do not appear in the diagonal matrix elements, so will 


not be adjusted. 


The NMREN/NMRIT method of calculation also leads 
to considerable difficulty in calculating meaningful 
error estimates for the final parameters. Ferguson and 


18b use a more elegant procedure than Swalen 


Marquardt 
and Nea but both calculate only the contribution 


of the errors in Fob to the errors in the parameters. 


Ss 
They ignore the additional errors from the residual 
differences between Hes obtained from (E.2), and the 
diagonal matrix elements of yD Sai anes with the 
final parameters. In principle, this error contribution 
could be included, but there is no obvious way to estimate 


the errors introduced by the deliberate exclusion of 


off-diagonal matrix elements from NMRIT. 
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All of these difficulties are avoided in the LAOCOON 
approach introduced by Castellano and Bothher=By, ~/4 which 
focusses attention directly on the relationship between 
transition frequencies and parameters. In the trial 
spectrum used for the assignment step, each calculated 
transition is accompanied by a single origin number which 
unambiguously designates the energy levels connected 
by the transition. However, when the analyst assigns 
an observed transition frequency to an origin number, 
he need not know which energy levels the transition connects. 
There is no need to include sufficient transitions to 
produce a fully connected energy level diagram, because 


there is no step analogous to NMREN. 


In the iterative step of LAOCOON, 17a 4 matrix 
Om pantwal, derivatives; 8V</8P; of the transition 
frequencies v, with respect to the NMR parameters p 
(chemical shifts and coupling constants), is calculated. 
These partial derivatives are obtained from the differences 
between appropriate eigenvalues of an approximate Hamiltonian 
(scemEq mlb ol jial® 


H 
approx 


= E.3 
av, /2P; dE / ops 9B /8Py ( ) 


The partial derivatives of the eigenvalues have been 


showns. @ to be identical to the diagonal elements of 


U* (Ha yprox/ 2P4)U: Thus partial derivatives of both 
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the diagonal and off-diagonal elements of H are 
approx 
calculated in the basis representation, and transformed 


using the eigenvectors. 


These partial derivatives are used to calculate 
corrections Ap to all of the parameters p being adjusted 

) FN NAS aie Ga Rae aed Hoe, 
The least squares solution to these equations gives the 


corrections to the parameters for use in calculating 


H approx for the next iteration. 


The LAOCOON method has been adapted for calculations 


17b 


using magnetic equivalence, and is) quite: suitable: for 


a program, such as NUMAR, that uses full factoring. | 
The transformation of the partial derivatives from the 
basis representation is readily done one matrix, at a 
time, and the matrix of coefficients which must be inverted 


to obtain the Ap values is only as large as the number 


of parameters. 

Because the LAOCOON least squares solution involves 
the entire Hamiltonian, and relates parameters to transitions 
directly, 1c produces “far more reliable estimates ior the 
errors in the parameters than does NMREN/NMRIT. 


r C.W. Haigh (private communication) has adapted the 


LAOCOON method of iteration to a program LACX, which 
uses X factoring and limited twofold frame symmetry. 


93 


ee 


ory 


SrOHpen 


¢ bewypplp 2 nie 
yt i oo i Lad & su aad apes 
eer Kom ia} eas 


wynd ST ea om Ms 


. f Raed 
, a. 
ig zaottawpe, saat /ad no lsubon aH 


“| 
ou oe stuctets, Yon! x0 rg th 


od fe: smote od, ot 


na ¥ 


3s) ‘a T bishti Th soussevire i wy * 
pried q eretsoneteq sina to nis we 


met 


r ( S16 nad ‘ 2 as" 


¥ 7 ® =~ 
¢ Ae 


‘ a) 
a. a. 
ved P 
ees Ci ie pee i 


7 Te ins 7 7 / 
j) Oa 
me esol 


_ 


A 


' ut a! wees. 
mia5 fh r gee ¥ w ara samet ee ) bitten ate and 2s ‘ 


to tented ren 


? : - : fi ‘an 
197 botrbhs- ceed ean badgem sab 
‘ hy 
. ? Fi unl | 
odiap 2i hire eos oovts {gv iupo eae 
. . E : a a 
waeT I : 2940 sate SAMUW ee. , do oe Sai 


ftravirso® D622" ‘er ots: ‘ge ‘noes BH rd 

‘ ogi 
Prenat orrd- sob. ibast mY big 20744 
» - % : et ‘ A : 
isitw ape oat 1003 Fo 3 cil bata , 


ri oS 


' ot ae aay 
, reer ny Lila ‘i : 
witulae eatempe panes | | 
i -pwtemsvor ey) heat i ghios Lime 
sthien tea “OLUREISS “TOK F wi i? 2hID 
| Aja esi pape iss ae 
y ae 
= eae tyes % r sean — 


= 
- 


va pay isa re 


| ‘iain 


i] 


sein te 


i ° 


i 
6 Vay le 


An iterative parameter adjustment step following 
the LAOCOON approach to iteration could be added to the 
NUMAR program with a minimum of difficulty. Most of the 
data required to calculate the parameter corrections 
(basis functions, eigenvalues, and eigenvectors) is already 
placed in auxilliary storage by the present program; 
observed transition frequencies, in order of line number, 
could. readily berplaced in auxilliary steragesin,the 


INPUT step. 


One could therefore add the parameter adjustment 
step as a completely independent subroutine (ADJUST in 
Fig. 4-1) with negligible changes in existing subroutines. 
ADJUST would require main memory space for a matrix of 
coefficients as well as for eigenvector matrices for 
a transformation similar to that in TRANS. It would 
therefore become the step limiting problem size for an 
iterative calculation; one could, if necessary, allow 


larger matrices, for non-iterative calculations. 
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APPENDIX F. MODIFICATIONS TO VARIAN SPECTROMETERS; 


INDOR METHODS. 


A Varian HA-60I spectrometer has been modified 
to permit INDOR spectroscopy 16and to make large lock 
offsets more convenient. Varian HA-100-12 and HA-100-15 
spectrometers have been slightly modified to permit more 
convenient use of large lock loffsetst All modifications 
are in the V4354 Internal Reference Stabilized NMR 
Controllers, which provide audio frequency magnetic field 
modulation, and detect signals at the modulation frequencies. 
They are, therefore, .effectivemfoneali radiofrequencies 


at which the spectrometers operate. 


Mipures!-U(a) Shows Saparciaieci rout diagram 
of the V4354 controller before modification. The Manual 
Oscillator can be set by hand to any frequency in the 
Gange 1500-3500 Hz., and the Sweep Oscillator, whose 
frequency is linked to recorder position, can sweep any 
of several frequency ranges all beginning at 2500 Hz., 
the widest extending to 3500 Hz. Sidebands of the basic 
radiofrequency (eg. 56.4 MHz. for '°F) are obtained by 
modulating the magnetic field at the audiofrequencies 
produced by these two oscillators. In Frequency Sweep 
mode, the Manual Oscillator sideband corresponds to the 
resonance frequency of the reference compound, and the 


Sweep Oscillator sideband is swept through the frequency 


am} LLaRORT eT 2 i. AEM, _s fz a1 nh PEAD 
aye mais | 


2 -oOt hy} bane, S . 


4 ce. ro 
xsrom SIS OF 


SEDI pari i. ‘fst pbon alt 23 By af anaes iaased ‘bam, nok 


7 « - 


nae -*: 


“ the 


x i 
dat & ae 
», =e 


‘ 


fy St #ehon stood vail 1414998 
£508 ‘Sart see or. a 


ree a flog obi aati, .eteftor 
$s 


cn ro 4 TROd Bact ae et oni. syste, swsobatods Mee 


* Lodeciigy arTono73 9898 ods de 


meveshh tio tec renga o: Cree (eyes oq 


‘ ¥ 
Cayo ott ~itod TH oak a noted: aeons deny. 
art agentes rms oO baad ta 2 of mo a 
saailw ~t62 ath 220 ah ot bos at 0828-003 
Ys. qoaw2 a 4 the i ‘sieo@ vabeiynen oS. ‘Bounit ang se 
sh 0025 te QaranZ aed its eager ‘eonsispar? tarew 
“ Oy , - I cn 
stend sid to rebasdebsae lt one, ‘oF < ee tee 
+ } . Re i | AS P 
yl Kies do ts (yet 0d | 


2otD tacipee tic LOB oils ts. 
weave. Konouperd tl 

pis ot gOapagoeres) ARAY 

TE DR, Pawan. am a 


yoaauper? oft Aguoriit: Apaas: € 


. Ay 7 2 a 
r : ry yi ; 7 7 a 7 
* . | 4 - - D f . _ vf . 
ia: ina i 


SWEEP MANUAL 
OSCILLATOR OSCILLATOR 


(REAR PANEL) 


MAN. OSC. 
ba LOCK 
‘PHASE CHANNEL 

REFERENCE 
SWEEP OSC. y SIGNALS SPECTRUM 
ouT ? ; 
z. (FRONT PANEL) CHANNEL 
MODULATION 


LEVEL CONTROLS 


ADDER AMPLIFIER 


(a) ORIGINAL CIRCUIT 
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(b) MODIFIED CIRCUIT 


FIG. F-l, MODIFICATIONS TO VARIAN V-4354 CONTROLLER. 
(SIMPLIFIED PARTIAL CIRCUIT DIAGRAMS) 


£4 


ad 
ashton 


{games png) 
“POT ae, eG 
Fe sORTHeD: abel 


Pat asi geen 


ol 
assis 
l Fae uiAe 


A ee 


Ca 


cael ray gt 
< ; 


ABIIOATHOD: a 
ic cates 2 


region whose spectrum is of interest. 


An audiofrequency phase shifter was inserted in 
the reference phase path of the spectrum phase detector, 
in both the HA-100-12 and HA-100-15 spectrometers. To 
offset the signal loss from the phase shifter the gain 
of the phase detector was increased. Otherwise, no wiring 
changes were made. When large lock offsets were used, 
the radiofrequency phase was adjusted for maximum lock 
Signal, and the audiofrequency phase was adjusted to 


give an absorption mode spectrum. 


Porethe spectra sat SA l@and@ lovey Nice) wenaprerr sy, 
tiesticernal manual"osei Plator-of "eaen -spectrome ter was 
dpsabled by “tire “removal “of "its “Oscrllator “Amplrtier*circuit 
card. The output of a Hewlett-Packard 5100B/5110B frequency 
synthesizer was fed in through the Man. Osc. Out jack 
Greco eet fay =) “CO=provide stapre, locker rrequencres up tu 
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The more thoroughly moditied circuit of-Fig..F=1(b), 
used in the HA-60I, allows a wider variety of operating 
modes. An external oscillator for large lock offsets 
is simply connected to Man. Osc. In; the phase shifter 
(as described above) allows correct adjustment of both 
lock and observing phases. For normal operation, essentially 
as in Fig. F-l(a), Man. Osc. Out is connected to Man. 

Osc. In, and Sweep Osc. Out to Sweep Osc. In. The third 


modulation input simplifies low power double resonance 


mt borveant aaw sertrde ia 
er rede saedg att bogs otis des 
ot | See 21- OE “nb 7 005. 
kas std. - 3S ont ite geena aay mot eg angie - 


- vas fx & 
aciréw or 52 Liecnadet bannstoid 286 wen 92. be: 


jaan aoe ob sethe tool sgret andi ae 


save 
be bag 
foo i ms ROR 70% bodes phe, anw zed 


_ 


7 eS = ' ae * «arr Saas ‘oboant simans 


(2 ,ataqnt) ish 0. @e2 bin ot. ne +8 asz0q2 ody 
may toodre Ane, te sodattéons! tauabni ‘Tan 
fieyrin wad Hi Gynt rote tiseO! 2$t Yo tevoner odt coe 

Ti 5 4 ai ‘PAHOOE. bretowT- trstwok-a. to rugine 8 
Lb s® 490 ‘woe cw wie aQuemese: mt bak paw 


sotainuper® Meet eidare: eee #. i 


; = a Ls * ll 


‘ a aes: oe 


, , ae ae 
ezeetto font raee 102, ee oe, fs rT = 


tottsd2 sated AF Pe yaa" ad Lesttt 329 fain | 
iged }0 Lyodvran tiie s297t0> ee ‘be 
re br eeeaS mek tetege earsote 3 | 
aft OF ‘bot Aithes ze 00 


ba isit ad Bil: 
Tak 


o ca 


a 

,Cdyf-4 .git tos kaoeia heer tboar cCaigusoat o10m ¢ AT 
. ] nani a | Pere 

ati tetAgo, to Wottey wobin 8 ewohts’ (108-4 of anid 


wo) 
9) 


pk 


oY hedemihe' aes Diet (oasupertotbhes ‘eit ‘b bat « 


: 15 
experiments. The use of the modified spectrometer for 


INDOR! is described below. 


For ally 56.4iMHzey) ?Byspéctrardé5.49 a Hewlett? 
Packard 5100B/5110B frequency synthesizer was connected 


to Man. Osc. In to provide stable lock frequencies up 
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For: the’ 56. 4° MHz." INDOR “spectra of §5.1,° the frequency 
Symenesizer~ was Set’ to 10640700"Hz) == With; the internal 
sweep’ Oscillator ‘connected’ to*Sweep"0sc.-iIn, a conventional 
spectrum of the =C-F region was obtained. The recorder 
wassaajusted to bring the-lowsfieid member of the central 
doublet into resonance, and the sweep oscillator frequency 
was measured using a Hewlett-Packard 5245L counter in 
period average mode. .The internal Manual Oscillator 
was then adjusted to this same frequency (about 2540 Hz.), 
and connected to Sweep Osc. In. It was discovered that 
ties Manual Oscillator hadia short "term stabvlity-of 20°02 Hz° 
(measured at 5 second intervals) and would typically 
gristsonly about.0./05°Hz.%in 1000 seconds. ©The modulation 
level for the observing channel was set to produce slight 


saturation. 


The INDOR sweep frequency v, was obtained by 
connecting a Data Royal F230B voltage controlled oscillator 
to the Double Resonance In jack. The oscillator was 
swept from 6800 to 6900 Hz. (i.e. through the -CF, 


resonances) in about 1000 seconds by its internal ramp 
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generator, while the Varian recorder was swept across 


the chart paper on its 1000 second sweep time setting. 


For 'H homonuclear INDOR experiments no external 
INDOR sweep oscillator is needed. With a suitable choice 
of lock frequency, the INDOR sweep range can be made 
to coincide with the sweep range of the internal sweep 
WoecurlelOr. Keter Inltidl setupmesimilar £0 CuaueaDoVe, 
Sweep Osc. Out is connected to Double Resonance In and 
Man. USc. Out to Sweep Osc. In for “the “actual “INDOR 


spectrum. 
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APPENDIX G CALIBRATION OF SPECTRA 


A basically similar procedure was followed in the 
calibration of all NMR spectra used in 85.1 and §5.2. 
At a number of calibration points across each spectrum, 
pen position was read to +0.01 cm. from the chart paper 
grid, and sweep frequency was read to +0.01 Hz. using 
a Hewlett-Packard 5245L counter in period average mode. 
These position/frequency pairs were used as input to a 
computer program which obtained a least-squares quadratic 
equation for frequency .as.a,function of chart position. 
Tne trequenciessof the spectral Llinesywere then obtained 


Brome their Chart positions, susine the calibration equation, 


Higher order polynomials were tried for some spectra, 
but in every case the, estimated. errors in the cubic and 
higher coefficients were much larger than the coefficients 
themselves. Fewer than five calibration points were 
avablablesfor some «spectra; »imitheseccases jhaslinear 


calibration equation was used. 


For the '°F single resonance spectra of §5.1, 
the calibration points were obtained by stopping the 
sweep between peaks while the spectrum was being recorded. 
The RMS deviation of the calibration points from the 


least squares curve was typically 0.05 Hz. 


For the ?°F INDOR spectra at 56.4 MHz., the sweep 


could not be stopped during frequency measurements. 
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During that time when a frequency measurement was being 
taken, the recorder pen was lifted from the paper to mark 
the location. The resulting calibration was estimated 

to be accurate to +0.1 Hz. Because of the poor signal 

to noise ratio in the INDOR spectra, the positions of 

the INDOR peaks were estimated to be accurate only to 
+0.15 cm. (+0.3 Hz.); averages from several INDOR spectra 
were used to identify corresponding single resonance 


peaks whose frequencies were better known. 


The methine proton spectra of the 2,5-dimethyl- 
3-hexenes were obtained without interrupting the sweep. 
A number of strong, well resolved peaks were then chosen 
as calibration points. Their frequencies were determined 
individually to 0.01 Hz., with the recorder stationary 
at the position of maximum pen deflection. Special care 
was taken to avoid saturation. At least two determinations, 
in which the peak was approached from opposite sides, 
were averaged for each calibration peak. This procedure 
nearly eliminated systematic errors resulting from the 
finite sweep rate in the original spectrum. The RMS 
deviation of calibration peaks from the least squares 


curve was typically +0.03 Hz. 


In the methyl and olefin regions of these molecules 
the frequency of every strong line was obtained directly, 
following the same procedure used for the methine calibration 


peaks. 
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